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Abstract: We report on the calculation of the next-to-leading order QCD corrections to the 
production of W-boson pairs in association with a hard jet at the Tevatron and the LHC, which 
is an important source of background for Higgs and new-physics searches. Leptonic decays of 
the W bosons are included by applying an improved version of the narrow-width approxima- 
tion that treats the W bosons as on-shell particles, but keeps the information on the W spin. 
Contributions from external bottom quarks are neglected, because they are either numerically 
suppressed or should be attributed to different processes such as Wt or tt production. A survey 
of differential NLO QCD cross sections is provided both for the LHC and the Tevatron. The 
QCD corrections stabilize the leading-order prediction for the cross section with respect to scale 
variations. However, the scale dependence of the next-to-leading order results for the LHC is 
only reduced considerably if a veto against the emission of a second hard jet is applied. In 
general, the corrections do not simply rescale the differential leading-order cross sections. In 
particular, their shapes are distorted if an additional energy scale is involved. 
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1 Introduction 



The search for new-physics particles — including the Standard Model Higgs boson — will be 
the primary task in high-energy physics in the era of the LHC. The extremely complicated 
hadron collider environment does not only require sufficiently precise predictions for new- 
physics signals, but also for many complicated background reactions that cannot entirely be 
measured from data. Among such background processes, several involve three, four, or even 
more particles in the final state, rendering the necessary next-to-leading-order (NLO) calcula- 
tions in QCD very complicated. This problem lead to the creation of an "experimenters' wishlist 
for NLO calculations" [1-3] that were still missing at that time, but are required for successful 
LHC analyses. The process pp — » W + W~+jet + X made it to the top of this list. Meanwhile 
the 2 — > 3 particle processes on the list have been evaluated at NLO QCD, and we further con- 
tribute to and extend that part in this paper; most notably also some of the 2 — > 4 processes 
have been calculated to NLO recently, viz. for the production of ttbb [4] and W + 3jets final 
states [5]. Moreover, benchmark results for the virtual corrections have been presented for a 
specific phase-space point for all 2 — » 4 processes on the list in Ref. [6]. 

The process of WW+jet production is an important source for background to the production 
of a Higgs boson that subsequently decays into a W-boson pair, where additional jet activity 
might arise from the production. In particular, it has been shown recently that the sensitivity 
of the Tevatron experiments CDF and DO concerning the discovery of a Standard Model Higgs 
boson in the mass range 135-190GeV could improve significantly by studying H — » W + W~ — » 
1 + 1 ~i>T together with an additional jet [7]. 

WW+jet production delivers also potential background to new-physics searches, such as su- 
persymmetric particles, because of leptons and missing transverse momentum from the W de- 
cays. Besides the process is interesting in its own right, since W-pair production processes 
enable a direct analysis of the non-abelian gauge-boson self-interactions, and a large fraction 
of W pairs will show up with additional jet activity at the LHC. Last but not least WW+jet at 
NLO also delivers the real-virtual contributions to the next-to-next-to-leading-order (NNLO) 
calculation of W-pair production, for which further building blocks are presented in Ref. [8]. 

Here we report on the calculation of the process pp/pp —* W + W~+jet + X in NLO QCD 
including leptonic W-boson decays. Results of this calculation on integrated cross sections, 
which are not sensitive to the W decays, have already been published in Ref. [9]. Parallel 
to our work, another NLO study [10] of pp — » W + W~+jet + X at the LHC appeared, where 
also the W decays have been taken into account. Moreover, a third calculation focusing on 
WW+jet production is still in progress [11]. In Ref. [3] the three different approaches are briefly 
described, and a detailed comparison of the virtual corrections at a single phase-space point is 
presented, revealing mutual agreement. 

The paper is organized as follows. In Section [2] we describe the general setup of the cal- 
culation, paying particular attention to the treatment of the various quark flavours. Section |3] 
provides details of the NLO calculation, and Section |4] describes the methods to include the 
W decays. Our detailed discussion of numerical results is presented in Section[5l which covers 
both integrated and differential cross sections. In Section[6]we give our conclusions, and the 
Appendix contains a derivation of the improved narrow-width approximation. 



1 




Figure 1: Diagrams contributing to one specific LO subprocess: The two W bosons couple 
directly to the fermion chain (upper diagrams) or by means of an intermediate neutral gauge 
boson (lower diagrams). 

2 Leading-order calculation and calculational framework 

At leading order (LO), hadronic WW+jet production receives contributions from the par- 
tonic processes qq — > W + W~g, qg — > W + W~q, and gq — » W + W~q, where q stands for up- or 
down-type quarks. Note that the amplitudes for q = u,d are not the same, even for vanishing 
light-quark masses. All three channels are related by crossing symmetry. The LO diagrams for 
the sample process uu — > W + W"g are shown in Figured! 



2.1 Quark-mixing effects 

To very good approximation the Cabibbo-Kobayashi-Maskawa (CKM) matrix can be as- 
sumed to be of block-diagonal form allowing mixing only between the two light generations, 
namely 





( V ud Vus 


()\ 




cos 6q sin 9q 0\ 


^CKM = 


V c d Vcs 







-sin^c cos 8q 







V 


V 


1/ 



where 6c denotes the Cabibbo angle. The approximation is justified by the fact that the ne- 
glected off-diagonal matrix elements are very small. Moreover, ignoring the masses of the light 
quarks, the dependence on the CKM matrix drops out for a remarkable set of subprocesses. For 
PP/PP ~* WW -i-jet + X, this happens if both W bosons couple to the same fermion chain, as 
illustrated in Figure [H Independent of gluonic couplings to this fermion chain, which do not 
affect the electroweak structure, the unitarity of the CKM matrix leads to 
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Figure 2: Electroweak part of diagrams with two W bosons coupling directly to the same 
fermion chain (left and central diagram): With light-quark masses neglected, all contributions 
with U U' or D ^ D' vanish in the block-diagonal approximation of the CKM matrix due 
to cancellations between diagrams with different intermediate states. The contributions with 
U = U' or D = D' behave as in the case of a trivial (unit) CKM matrix. Diagrams with an inter- 
mediate neutral electroweak vector boson (right diagram) do not depend on the CKM matrix, 
and thus Q = Q' holds. 



Figure 3: A sample diagram with two W bosons coupling to different fermion chains: Contri- 
butions of this kind yield a non-trivial dependence on the CKM matrix elements. 

with the nomenclature of Figure [2l when the intermediate quark state is summed over. Diagrams 
with the W-boson pair coupling to the fermion chain by means of an intermediate vector boson 
are independent of the CKM matrix anyway. Therefore, a remarkable set of subprocesses of 
pp/pp — » WW+jet + X is not influenced by the explicit entries of the CKM matrix due to its 
unitarity. 

The only subprocesses contributing to pp/pp — » WW+jet + X that depend on the explicit 
entries of the CKM matrix are those containing two fermion chains with the two W bosons 
coupling to different fermion chains at least in some diagrams. These are the real-emission 
subprocesses including both two external up-type and two down-type (anti-)quarks. A sample 
diagram is shown in Figure [3] 

2.2 Treatment of bottom (anti-) quarks 

The situation of a hard jet resulting from an outgoing bottom (anti-)quark has to be consid- 
ered with care: The long bottom lifetime results in a resolvable second vertex that is displaced 
from the primary interaction point, allowing in principle to tag the bottom flavour. This, how- 
ever, does not imply that all events involving outgoing bottom quarks can be separated. For 
instance, in real-emission subprocesses only one of the two outgoing jets must be detected: 
This can be the light jet as well, while the b-jet leaves the detector unseen in direction of the 
beam axis. Moreover, the b-tagging efficiency is below 100%, so that not all "b-jet events" 
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can be isolated. The influence of incoming bottom (anti-)quarks on the hadronic cross section 
is suppressed with respect to other incoming (anti-)quarks by the small bottom PDFs in the 
colliding hadrons. 

In the process class pp/pp — » WW + jet + X, however, the suppression of subprocesses with 
external bottom (anti-)quarks may be overcompensated by top resonances, always showing up 
along with final- state bottom (anti-)quarks. Moreover, outgoing bb pairs appearing in real- 
emission subprocesses even contain two potentially resonant top-quark propagators in some 
diagrams without PDF suppression. These subprocesses should, however, in general not be 
assigned to WW+jet production, since they are actually off-shell continuations of tt production. 
From this point of view, the subprocesses with only one outgoing bottom (anti-)quark can be 
seen as contributions to W"t production or W + t production, with the off-shell decay of the 
top (anti-)quark included. In our approach, precisely described in the following, we exclude 
these resonance reactions upon omitting contributions from channels involving external bottom 
(anti-)quarks. Thus, if desired, contributions from these channels could be easily added, but 
their calculation requires a consistent treatment of top (anti-)quarks as unstable particles, as for 
instance provided by the complex-mass scheme [19, 20]. 

For the calculation of cross sections, we follow two different strategies regarding the treat- 
ment of bottom (anti-)quarks: 

( i) Five-flavour scheme 

In the first approach, bottom (anti-)quarks are treated as massless particles, so five-flavour 
PDFs and a five-flavour running of a s are used. In this framework, the strategy is to neglect 
all contributions containing external bottom (anti-)quarks. For the initial state, this approach 
is justified by the smallness of their PDFs. Final-state bottom (anti-)quarks are excluded by 
the assumption that their signal can be distinguished from that of a light-quark jet by means of 
b-tagging. 

An advantage of this procedure can be seen in the fact that the influence of the off-shell 
continuations of tt, W"t, and W + t production, which is explained in the foregoing passage, is 
simply left away. The respective contributions should be added to dedicated calculations for 
these processes. Note, however, that the naive application of this procedure would lead to an 
ill-defined cross section containing mass singularities. The problem is due to the g — » bb split- 
ting which by anti-b-tagging would be removed from the real corrections. This contribution is 
required to cancel the mass singularities from the bottom-quark loop in the gluon self-energy 
contributing to the virtual corrections. We treat this case similar to what has been done in 
Ref. [21] and is also discussed in Ref. [22]. If the two bottom quarks are combined according 
to the applied jet algorithm the resulting jet carries no (net) b-charge and is thus counted as a 
"light" jet. On the other hand if the two bottom quarks form two individual jets we assume 
them as tagged and ignore this configuration. Applying this procedure restricts the phase-space 
integration to the collinear configurations which are needed to cancel the mass singularities in 
the virtual corrections. Within the dipole subtraction formalism we combine the corresponding 
integrated dipoles with the virtual corrections to obtain a finite result. In addition one obtains 
a contribution from the real matrix element with the unintegrated dipoles subtracted. However 
this contribution is only integrated over the collinear phase space region. Since in this region 
the unintegrated dipoles approximate the real matrix elements we can ignore this contribution 
to very good approximation. The case in which the bottom-quark pair is not produced via gluon 
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splitting is highly suppressed by anti-b-tagging. The neglected contribution only delivers size- 
able contributions from diagrams including resonant top quarks. Exactly these contributions, 
however, are meant to be left out, because they should actually be assigned to different process 
classes. 

( ii) "Four-flavour scheme " 

In the second approach, the bottom (anti-)quarks are understood as massive particles. We 
use four-flavour PDF's to describe the parton contents of the proton. Here, no top resonances 
show up, since no bottom-(anti-)quark densities are taken into account and no mixing between 
the two light and the third generation takes place in the chosen approximation of the CKM 
matrix. As a consequence no single outgoing bottom (anti-)quarks appear. No bb pairs are 
taken into account in final states which is justified by assuming anti-b-tagging. The running 
of a s is driven only by the four remaining light quarks with both bottom- and top-quark loops 
in the gluon self-energies subtracted at zero momentum. As explained in Section [31 no large 
corrections arise from terms proportional to a s lnmb, because counterterms always contribute 
in the combination j6Za + 8Z gs . Thus, the a s lnmb term from the renormalization of the strong 
coupling cancels against the corresponding term from the wave-function renormalization of the 
gluon. In this scheme, such a cancellation always takes place if the number of external gluons 
is equal to the number of strong couplings in the considered LO process. 

3 Details of the NLO calculation 

In order to prove the correctness of our results we have evaluated each ingredient twice using 
independent calculations based — as far as possible — on different methods, yielding results in 
mutual agreement. 

3.1 Virtual corrections 

The virtual corrections modify the partonic processes that are already present at LO. At 
NLO these corrections are induced by self-energy, vertex, box (4-point), and pentagon (5-point) 
corrections. For illustration the pentagon graphs, which are the most complicated diagrams, are 
shown in Figure |4] for one partonic channel. At one-loop level WW+jet production also serves 
as an off-shell continuation of the loop-induced process of Higgs+jet production with the Higgs 
boson decaying into a W-boson pair. In these diagrams the off-shell Higgs boson is coupled 
via a heavy-quark loop to two gluons; the graphs for this mechanism are shown in Figure [5] 
together with vertex-correction diagrams with an intermediate electroweak vector boson. The 
Higgs resonance is included in our calculation upon employing a fixed decay width Th, i.e. by 
replacing the propagator denominator (p^ - M^) by (pj-j - + iMnYyi), where pn and Mh 
denote the momentum and the mass of the Higgs boson. The box contributions to the fermionic 
corrections are depicted in Figured 

An on-shell renormalization is performed for the wave functions of the external QCD par- 
tons and an MS renormalization for the strong coupling with the massive-quark loops in the 
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Figure 4: Pentagon diagrams for uu — > W + W g. 
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Figure 6: Fennionic box diagrams for qq — » W + W g. 



gluon self-energy subtracted at zero momentum. Since only massless external quarks are in- 
volved, the whole counterterm amplitude can be written as 



Met = (sZq + l -8Z A + 8Zg\ M LO , 



where the renormalization constants are calculated as follows, 
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The sum over q runs over all massive quarks, namely q = t in the five-flavour scheme (TVf = 5) 
and q = b, t in the four- flavour scheme (JVf = 4). Both UV and IR (soft and collinear) divergences 
are regularized in D = 4 - 2e dimensions, and the divergences arising as poles in e are quantified 
by 

'% 2 V r d+e) 



M 2 W 
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(3.5) 



where n is the arbitrary reference scale of dimensional regularization. Superscripts "UV" or 
"IR" on A;t(/i) indicate the origin of the divergences. 

As is well known, the use of the Dirac matrix 75 deserves some care within dimensional 
regularization, because the Dirac algebra with 75 does not admit a straightforward analytical 
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continuation to D ± 4 dimensions. In both calculations of the virtual corrections we use a simple 
recipe with a 75 that anticommutes with all other Dirac matrices but nevertheless obeys the usual 
trace relations from four dimensions [j] Although this approach is, of course, not fully consistent 
from the mathematical point of view, it is well known that it delivers the correct results at NLO 
as long as closed fermions loops are calculated "generically", i.e. they are calculated for general 
fermion flavours followed by a subsequent summation over all fermions by only changing the 
quantum numbers. The arguments supporting this scheme in the presence of UV divergences 
can be found in Ref. [41] (see also references therein). In view of dimensionally regularized IR 
divergences, the situation is much simpler because of their factorization off tree-level structures. 
As long as it is guaranteed that these residual tree structures are treated in the same way in the 
virtual and real corrections, between which the cancellations of IR divergences take place, the 
IR-finite sum is well defined. In fact for the compensation of IR divergences no (anti-)commu- 
tation rules involving 75 are needed at all. Since the naive scheme described above does not 
induce any spurios terms from the regularization, it automatically delivers the correct IR-finite 
sum of virtual and real corrections. 

Version 1 of the virtual corrections is essentially obtained as for the related processes of 
ttH [23] and tt+jet [24] production. Feynman diagrams and amplitudes are generated with 
FeynArts 1 .0 [25] and further processed with in-house Mathematica routines, which automati- 
cally create an output in Fortran. The IR divergences (soft and collinear) are analytically sepa- 
rated from the finite remainder in terms of triangle subdiagrams, as described in Refs. [23,26]. 
This separation, in particular, allows for a transparent evaluation of so-called rational terms 
that originate from Independent terms multiplying IR divergences, which appear as single or 
double poles in e. As generally shown in Ref. [4], after properly separating IR from UV di- 
vergences such rational terms originating from IR divergences completely cancel; this general 
result is confirmed in our explicit calculation. For the results presented in Ref. [9], the pen- 
tagon tensor integrals were directly reduced to box integrals following Ref. [27], while box and 
lower-point integrals were reduced a la Passarino-Veltman [28] to scalar integrals. This proce- 
dure completely avoids inverse Gram determinants of external momenta in the reduction step 
from 5-point to 4-point integrals, but the reduction of box and lower-point tensor integrals in- 
volves such inverse determinants via the Passarino-Veltman algorithm. Although these inverse 
determinants jeopardize the numerical stability in regions where such determinants are small, 
sufficient numerical stability was already achieved. Meanwhile the tensor reduction has been 
further improved using the methods of Ref. [32]. In detail the reduction of pentagons is per- 
formed by the more recent procedure of Ref. [32] (similar to a method proposed in Ref. [33]), 
and the Passarino-Veltman reduction for 4-point integrals, etc., is supplemented by the dedi- 
cated expansions for small Gram and kinematical determinants in the regions where these de- 
terminants become smallH The scalar one-loop integrals are either calculated analytically or 
using the results of Refs. [29-31]. 

Version 2 of the evaluation of loop diagrams starts with the generation of diagrams and 
amplitudes via FeynArts 3.4 [35] which are then further manipulated with FormCalc 6.0 [36] 
and eventually automatically translated into Fortran code. The whole reduction of tensor to 

The authors of Ref. [11] adopted the more rigorous, but cumbersome approach of 't Hooft and Veltman [40], 
where 75 is split into an anticommuting four-dimensional part and a commuting remainder, and find full agreement 
with our result (see also Ref. [3] for some details). 

2 Similar procedures based on expansions in small determinants have also been proposed in Ref. [34]. 
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scalar integrals is done with the help of the LoopTools library [36], which employs the method 
of Ref. [27] for the 5-point tensor integrals, Passarino-Veltman [28] reduction for the lower- 
point tensors, and the FF package [37, 38] for the evaluation of regular scalar integrals. The 
dimensionally regularized soft or collinear singular 3- and 4-point integrals had to be added 
to this library. To this end, the explicit results of Ref. [26] for the vertex and of Ref. [39] for 
the box integrals (with appropriate analytical continuations) are taken. Actually the FormCalc 
package assumes a four-dimensional regularization scheme for IR divergences (such as the 
concept of an infinitesimal photon mass in QED), i.e. rational terms of IR origin are neglected 
by FormCalc. However, as mentioned above, in Ref. [4] it was generally shown that such 
rational terms consistently cancel if UV and IR divergences are properly separated. Owing to 
this property the algebraic result of FormCalc for the unrenormalized amplitudes could be used 
without any modification, apart from supplementing the needed IR-singular scalar integrals. 

3.2 Real corrections 

The matrix elements for the real corrections are given by the processes — » W + W _ qqgg and 
—* W + W~qqq'q' with a large variety of flavour insertions for the light quarks q and q'. The 
partonic processes are obtained from these matrix elements by all possible crossings of quarks 
and gluons into the initial state. The evaluation of the real-emission amplitudes is performed in 
two independent ways. Both evaluations employ (independent implementations of) the dipole 
subtraction formalism [42] for the extraction of IR singularities and for their combination with 
the virtual corrections. 

Version 1 employs the Weyl-van-der-Waerden formalism (as described in Ref. [43]) for the 
calculation of the helicity amplitudes. The phase-space integration is performed by a multi- 
channel Monte Carlo integrator [44] with weight optimization [45] written in C++, which is 
constructed similar to Racoon WW [ 1 9,46] . The results for cross sections with two resolved hard 
jets have been checked against results obtained with Whizard 1.50 [47] and Sherpa 1.0.8 [48]. 
Details on this part of the calculation can be found in Ref. [49]. In order to improve the integra- 
tion, additional channels are included for the integration of the difference of the real-emission 
matrix elements and the subtraction terms. 

Version 2 is based on scattering amplitudes calculated with Madgraph [50] generated code. 
The code has been modified to allow for a non-diagonal quark mixing matrix and the extraction 
of the required colour and spin structures. The latter enter the evaluation of the dipoles in the 
Catani-Seymour subtraction method. The evaluation of the individual dipoles was performed 
using a C++ library developed during the calculation of the NLO corrections for tt+jet [24]. For 
the phase-space integration a simple mapping has been used where the phase space is generated 
from a sequential splitting. 

In the Catani-Seymour dipole subtraction formalism all the IR divergent pieces with LO 
kinematics are collected in the so-called J-operator — for details we refer to Ref. [42]. Since 
the IR finiteness of the virtual corrections combined with the IR singularities obtained from the 
real corrections provides an important check of the calculation, we reproduce here the explicit 
form of the J-operator, 

I(p a ,Pb\---,Pi) = J^dcr^ b (p a ,p b ;...pi) = -^\^LO,ab(Pa,Pb\---,Pi)\ 2 
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where 3i\jo,ab denotes the colour- stripped LO amplitude for the incoming quark-antiquark pair 
ab and the outgoing gluon i. Due to the simple colour structure the J-operator does not involve 
non-trivial colour correlations. The J-operators for the remaining subprocesses are obtained by 
exchanging the respective momenta and using the crossed matrix elements. 



4 Inclusion of gauge-boson decays 



Since the produced W bosons are unstable particles, their decays should be included into 
the analysis. In this context, especially the case of both gauge bosons decaying leptonically 
is of interest due to its clean signature in the detector. Therefore, only leptonic gauge-boson 
decays are considered here. The inclusion of the decays is performed following three different 
strategies for the LO processes, which are presented in the following paragraphs. The improved 
narrow-width approximation — the on-shell approximation together with the exact treatment of 
the W-polarization — delivers an appropriate compromise between the complexity of the calcu- 
lations and the accuracy of the results, which is demonstrated in a comparison of LO results 
presented in Section l5.3.1l Therefore, the NLO calculation is performed by means of the im- 
proved narrow-width approximation. 

4.1 Full calculation with off-shell W bosons 

The diagrams contributing to pp/pp — > t~t'vi>vi + jet + X can be subdivided into two classes, 
namely diagrams showing two gauge-boson propagators that can become resonant and other 
diagrams containing only one W resonance. The doubly-resonant diagrams comprise all the 
WW+jet-production diagrams of Figure[7] with the two massive gauge bosons decaying into 
leptons. 

Aside from this doubly-resonant group, a number of diagrams that in general contain only 
one possibly resonant gauge-boson propagator contributes. Their topologies are characterized 
by the situation that one neutral gauge boson decays to a lepton-antilepton pair with the second 
pair attached to this leptonic fermion chain via a W boson. All diagrams for this configuration 
are given in Figure [8] for one specific partonic channel. The configuration with both the Z and 
the W propagator being resonant is kinematically allowed, but does not significantly contribute, 
since the corresponding phase-space region is too small. 

The singularities appearing at the poles of the LO gauge-boson propagators are treated by 
introducing decay widths in the propagators. More precisely, the gauge-boson widths are intro- 
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Figure 7: Diagrams with two resonant W propagators in the partonic subprocess 



uu 



duced in terms of complex gauge-boson masses according to the complex-mass scheme at LO, 
which is described in Refs. [19,20]. Thus, in the LO amplitudes we perform the substitutions 



M w -» M w - ir w M w , Ml^Ml- iY z M z . 



(4.1) 



To preserve gauge invariance, the gauge-boson masses have to be treated as complex quantities 
everywhere, in particular in the definition of the weak mixing angle, 

M w M w -ir w M w 
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which renders all couplings complex that are derived from this quantity. The values of the 
gauge-boson widths are calculated at NLO QCD level with vanishing fermion masses, 
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where /V c = 3 is the number of quark colours and gy/y/ denote the chiral couplings of V = 
y,Z,W to the fermions /,/', 
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Figure 8: Diagrams with only one resonant W propagator in the partonic subprocess 
uu -> VeeVv^g. 

Here Qf and = +j axe the relative electric charge and the third component of the weak 
isospin of /, respectively. In (14.31 ), the sums run over all three lepton generations and the 
two light-quark families, in (14.41) over the three charged leptons and neutrinos and the five light 
quarks. The W and Z decay widths are calculated using real masses Mw and Mz, and the strong 
coupling is set to the measured value a s (Mz) quoted below. Further improvements would go 
beyond NLO QCD accuracy. 

4.2 Narrow-width approximation 

In the naive narrow-width approximation (NWA), the produced W bosons are treated as 
on-shell particles in the production process. The leptonic decays are assumed to be isotropic 
in the rest frames of the respective W bosons, because their spin information from the produc- 
tion process is dropped. Therefore, the squared matrix elements of the WW+jet-production 
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Figure 9: NWA diagrams with two on-shell gauge-boson propagators for the partonic process 
uu — > v e e + /i - v^g at LO. The on-shell propagators are denoted by slashed propagator lines. 

subprocesses can be used without modifications by only multiplying them with the respective 
branching ratios for the W-boson decays, which are given by 

RP r W+-+vil + R „ r W-^l-vi 

±SK w +^ yi i+ - , ±SK W — — >1— vi - • (4.0) 

1 w 1 w 

The partial widths for leptonic W decays are 

rw+^vii+ = Iw--»i-y, = — Mw(gwvii) ' (4.7) 

since the leptonic decays do not receive NLO QCD corrections. 

Consequently, the approximation of on-shell gauge bosons restricts the diagrams contribut- 
ing to the subprocesses to those containing two resonant W propagators, which stand in one-to- 
one correspondence to the pure gauge-boson pair-production diagrams in Figure Ul The analo- 
gous NWA diagrams are collected in Figure [9] for one LO subprocess. 

4.3 Improved narrow-width approximation 

In the comparison that is presented in Section l5.3.1l for LO calculations, the naive NWA 
turns out to be not satisfactory. While the integrated cross sections show modest deviations 
between the full amplitude calculation and the naive NWA, large discrepancies arise in some 
phase-space regions if differential cross sections are considered. 

An appropriate compromise between the full calculation and the naive NWA is obtained 
by an improved NWA (iNWA). Here, the gauge bosons are still treated as on-shell particles, 
but their spin information is kept. It is used to improve the description of the leptonic W- 
boson decays, which are not isotropic in the respective W-boson rest frames owing to the V-A 
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structure of the decay. Details on how to keep the spin information of the decaying gauge 
bosons are provided in App. [Al 

The diagrams that are relevant for the iNWA are the same as for the NWA, which are shown 
in Figure [9] for one specific channel in LO. The modification of amplitudes is essentially the 
same for LO subprocesses and all contributions to the NLO QCD cross section: the polarization 
vector of each of the outgoing gauge bosons is replaced by the leptonic currents of its decay 
products. For the rest of the production amplitude, the momenta of the gauge bosons are set 
on shell. The Breit-Wigner propagators arising in the absolute square of the amplitudes are 
separated. They are replaced by delta functions with an appropriate normalization obtained by 
integrating over the Breit-Wigner propagator in the limit Ty — > 0. 

Since the iNWA turns out to reproduce the full calculation to sufficient accuracy in the LO 
comparison of Section 15". 3. 11 this strategy for describing the leptonic decays is applied for the 
NLO QCD calculations. 

5 Numerical results 

5.1 Setup and input parameters 

For the numerical evaluations, the following SM parameters [51] are used, 

M w = 80.425GeV, M z = 91.1876GeV , = 1.16637 x l(T 5 GeV _2 , 

M H = 150GeV, r H = 0.017GeV, m t = 174.3GeV. (5.1) 

The total Higgs-boson width Th has been calculated with Hdecay [52]. The only contribution 
from a Higgs boson at NLO QCD is given by the loop diagams in Figure [5J Since these diagrams 
turn out to contribute less than 0.3% to the cross sections for the chosen Higgs mass, their impact 
is not investigated further^] The electromagnetic coupling a is evaluated from (15.11 ) via 



V2G^ (l - 



Ml 



Mi , 

a = a G ^ = ^ ^- . (5.2) 

This choice of a absorbs some universal corrections to the electroweak coupling (running of 
a to the electroweak scale, leading corrections to the p-parameter), as is explained, e.g., in 
Ref. [53]. The widths of the weak gauge bosons are calculated according to Eqs. ( 14.31) and 
(14.41) — in case of the production of stable weak gauge bosons they are, of course, set to zero — 
from all decay channels at NLO QCD with the fermion masses neglected. With the value of the 
strong coupling at the scale Mz taken from Ref. [54], 

a s (M z ) = 0.1176, (5.3) 

the calculated widths are 

T w = 2.0996 GeV , Y z = 2.5097 GeV . (5.4) 



3 For Mh > 2Mw the included Higgs-boson effects account for the interference between a H(— > WW) + jet signal 
and the coherent irreducible WW + jet background. A proper description of the signal requires, of course, also the 
squared resonance diagrams and contributions from gg fusion with the corresponding radiative corrections. 
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The SM parameters already used in the publications on WW+jet production [3, 9, 55, 56] 
are not replaced by the more recent values [54] in order to facilitate comparisons; the numerical 
impact of the slightly changed measured value of Mw is negligible anyway at the required 
accuracy. 

The values for the strong coupling in the amplitude calculation are evaluated according to a 
1 -loop-running at LO and a 2-loop-running at NLO as described in Ref. [57], 

tfsj-loop^ren) = ~3 — > (5-5) 



12 * m A^c D 



°'s,2-loop(A'ren) — ( ^s,l-loop(A'ren) 



6(153 -19/Yf) V ^ 2 qcd) 
(33-2/V f ) 2 m j4n. 



(5.6) 



\ A QCD J 



where // ren is the renormalization scale, /Vf the number of active light-quark flavours, and Aqcd 
the QCD scale parameter. The values of Aqcd are chosen as prescribed by the applied PDF sets: 
In the five-flavour scheme, the PDFs of CTEQ6 [58,59] are used with Nf = 5, namely CTEQ6L1 
with Aqcd = 165 MeV at LO and CTEQ6M with A QC d = 226 MeV at NLO. In the four-flavour 
scheme, the PDFs of MRST2004 [60] with Nf = 4 are taken, namely MRST2004F4LO with 
Aqcd = 220 MeV at LO and MRS T2004F4NLO with Aqcd = 347 MeV at NLO. The renor- 
malization of a s is performed as described in Section [3TTT i.e. with the heavy-quark loops in the 
gluon self-energy decoupled: In the five-flavour scheme only the top-quark loop is decoupled, 
in the four-flavour scheme both the top- and the bottom-quark loops are decoupled. 
The Cabbibo angle in the approximation (12.11) for the CKM matrix is set to 

6 C = 0.227 , (5.7) 

the explicit entries are calculated from this. 

To give an IR-safe definition of the cross section we apply the successive-combination jet 
algorithm of Ref. [61] with R = 1 to decide whether two final- state partons can be resolved as 
two separated jets or whether they have to be combined to only one jet. A dependence on the 
specific jet algorithm only arises through the real-emission subprocesses, because the LO and 
all other NLO contributions contain only one parton in the final state which is identified with the 
jet. To resolve the additional hadronic jet, a cut on the transverse momentum is applied. This 
renders the cross section also IR safe which would otherwise diverge already at LO. Different 
values are used for /?T,jet,cut which are explicitly given for the respective results. 

If the weak-gauge-boson decays are included, a set of additional cuts is applied which is in 
general not necessary for the finiteness of cross sections, but provides results that are closer to 
the experimental situation. These additional cuts are 

|7/jetl < 4.5 , PTJepton > 25 GeV , |?7leptonl < 2.5 , (5.8) 

where rj is the pseudo-rapidity and pj the transverse momentum of the respective particle. All 
leptonic cuts are applied only to charged leptons, of course. As an additional condition the 
missing transverse momentum, which is deduced from the sum of the neutrino momenta, is 
required to obey 

p T ,m 1S s>25GeV. (5.9) 
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The missing transverse momentum is particularly of interest if WW+jet is considered as a back- 
ground process for SUSY searches. Here, the missing transverse momentum mimics the lightest 
supersymmetric particle that also leaves the detector undetected. 

Furthermore, isolation cuts are applied in order to separate all visible leptons from each 
other and from hadronic jets, 

^leptonjet > 0.4 , ^lepton.lepton > 0-2 , (5.10) 

where R a i, = a/ (A(p a b) 2 + (An a b) 2 . The angular distance between the two particles a,b in the 
transverse plane is denoted by A<p a b, and An a b is the difference of their pseudo-rapidities. 

5.2 Results on integrated cross sections 

We start our discussion of numerical results with integrated cross sections. In order to render 
these cross sections insensitive to the decay of the W bosons, we do not apply the cuts on the 
decay leptons specified in Eqs. (15.8I) -( 15.10I ) in this section. 

5.2.1 LO cross sections 

The LO cross sections get contributions from qq, qg, and gq initial states. In this section, 
the relevance of the respective partonic channels is discussed both for results on proton-proton 
collisions at yfs = 14TeV (LHC setup) and on proton-antiproton collisions at yfs = 1.96TeV 
(Tevatron setup). Scale variations of a factor 10 around the central scale, which is chosen to be 
Mw, are considered. Here and in the following sections, the common scale /i = // ren = yUfact is 
used, i.e. renormalization and factorization scales are set equal and are varied simultaneously. 

Figure [10] shows the results for the various partonic contributions to WW+jet production, 
evaluated in the five-flavour scheme. At the LHC, each partonic process involves at most one 
valence quark. Therefore, the gluon flux, which is essentially larger at the LHC compared to 
Tevatron, in general leads to larger contributions from qg channels compared to qq annihilation, 
whereas the gq initial states contribute significantly less since no valence partons are involved 
here. The situation at Tevatron is different, because the proton-antiproton collisions provide qq 
contributions with two valence partons that dominate the cross sections. The scale dependence 
of the cross sections at Tevatron turns out to be stronger than at LHC, which is due to the 
factorization- scale variation. While the renormalization- scale variation only reflects the running 
of the strong coupling a s , which is the same in both cases, the factorization scale dependence 
is quite flat for nearly all channels at LHC in the considered range. At Tevatron, however, an 
increase of roughly the same order as the one arising from the renormalization- scale dependence 
is found in the direction of lower scales. A detailed analysis of the scale dependence of LO cross 
sections to WW+jet production is provided in Ref. [49]. 

In case of the LHC, the partonic contributions with up-type quarks dominate over those 
with down-type quarks, which can be understood by comparing the valence-quark PDFs in the 
protons. The antiquark PDFs of the lightest generation result in an inverted order of the gq chan- 
nels at the LHC where the gd channel prevails the gu contribution. For Tevatron, the described 
valence-PDF effect occurs twice in the qq channel, leading roughly to a factor 4 between uu 
and dd. The channels qg and gq deliver exactly the same contribution to the integrated cross 
sections due to the charge-conjugation-invariant hadronic initial state. Naturally, the channels 
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Figure 10: Partonic contributions to the LO cross sections of WW+jet at the LHC (l.h.s.) and 
the Tevatron (r.h.s.): The variation of the common scale \x = jUf ac t = y^ren is shown for all chan- 
nels. If more than one curve with identical colours is given, the first one mentioned in the key 
corresponds to the upper curve. The contribution named "total" does not contain contributions 
from external bottom quarks. 



involving valence-(anti-)quarks dominate, and therefore the described PDF effect leads again to 
larger contributions from ug/gu compared to dg/gd. 

In addition, for both colliders the channels involving bottom flavours are shown in Figure [101 
As expected, bb annihilation is numerically negligible in both cases, whereas the bg channel 
and the gb channel — whose contributions to the integrated cross section are the same — each 
account roughly for the same amount as all the subprocesses involving (anti-)quarks of the 
two light generations at the LHC. As already explained in Section [2T2l this is due to the fact 
that these subprocesses actually describe resonant and non-resonant W~t and W + t production, 
respectively, with the top decays included. In addition the bottom PDF is due to the gluon 
splitting g — » bb. The bg-flux is thus rather large at the LHC due to the large gluon lumionsity. 
As mentioned earlier the corresponding contributions should be attributed to different processes, 
therefore these channels and the respective NLO corrections are not taken into account. For 
Tevatron, their numerical impact is not that large, since the centre-of-mass (CM) energy to 
produce top resonances in addition to a W boson is rarely available and the gluon flux is small. 
The respective channels are, however, treated in the same way as for the LHC. The different 
shapes of the scale-variation curves with bottom flavours result from a strong decrease of the 
bottom PDFs for smaller values of the factorization scale. 

5.2.2 NLO QCD cross sections 

After analyzing the individual contributions to the LO cross sections in the previous section, 
the effects of NLO QCD corrections on the WW+jet cross sections are discussed here. For 
the numerical results of this and the following sections, two different definitions of the NLO 
observables are used. The one observable is defined more inclusively by only requiring at least 
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A«/Mw /i/M w 

Figure 11: Comparison of WW+jet production cross sections in the LHC setup with 
PTjet > 50GeV and for Tevatron with /?T,jet > 20GeV: The straight lines show the results cal- 
culated with the five-flavour PDFs of CTEQ6, the dashed lines those calculated with the four- 
flavour PDFs of MRST2004F4. Contributions from external bottom (anti-)quarks are omitted, 
as described in Section [2T2l 



one hard jet with a minimum transverse momentum after application of the jet algorithm. The 
more exclusively defined observable applies a veto on a second separable hard jet and describes 
therefore genuine WW+jet production. To this end, real-correction events with two jets that 
fulfill the /?T,jet,cut condition and are not combined by the jet algorithm are not counted in the 
more exclusive observable. The phase-space regions in the real-correction subprocesses that are 
relevant for curing infrared singularities from the virtual corrections are not influenced, since the 
applied restriction only refers to genuine WW+2jets events. Therefore, the difference between 
the results for the two NLO observables is precisely given by the respective LO observable of 
WW+2jets production — evaluated, however, in the NLO setup (NLO PDFs, 2-loop running of 

( i) Four-flavour versus five-flavour scheme and scale dependence of NLO QCD cross sections 

We start the discussion of the scale dependence of the NLO cross sections with a compari- 
son between the results obtained in the four-flavour and the five-flavour schemes, respectively, 
as described in Section [2T2l The NLO QCD cross sections for WW+jet production are pre- 
sented in Figure QT] for the LHC setup with pTjet,cut = 50GeV and for the Tevatron setup with 
£>T,jet,cut = 20GeV, both for the four-flavour calculations with MRST2004F4 PDFs and the five- 
flavour calculations with CTEQ6 PDFs. In the four- flavour calculation, loop diagrams involv- 
ing (massive) 3rd-generation quarks are included. However, their contributions account for less 
than 0.2% to the integrated cross sections. 
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In case of the LHC, the relative deviations between the NLO cross sections in the two 
schemes behave as follows: while the deviation significantly increases in direction of lower 
scales at LO, which leads, however, to relative deviations of less than 3% in the whole range, 
it is nearly flat for both NLO observables where the results deviate by less than 2%. For the 
Tevatron, the comparison between the two approaches reveals deviations of less than 2% at LO 
over the full range of scale variations. The NLO cross sections turn out to be slightly larger 
in the CTEQ6 approach, but on a level of less than 4% for the whole range of interest. The 
smallest scales depicted may be ignored, because they obviously describe a region where the 
NLO calculation is not a good approximation anymore, as one can see from the steep decline 
of the absolute values for the NLO cross sections in that range of scales. 

Since the results evaluated with different PDF sets can in general not be expected to agree 
exactly — not only in the present case where different descriptions of the bottom (anti-)quark are 
applied — this analysis can be seen as a confirmation that the approximation that is applied in 
the five-flavour calculation delivers fully acceptable results. Moreover, the deviations between 
the two approaches are far below the expected experimental errors. The good agreement has 
been further confirmed by performing the same comparison for />T,jet,cut = 

lOOGeV at LHC 

and for /?T,jet,cut = 50GeV at Tevatron, yielding the same level of agreement. In addition, also 
differential cross sections have been compared with no significant deviations showing up in any 
phase-space regions. Numerical results of these further comparisons are omitted here, because 
they are not supposed to provide new insight. In the following, only the five-flavour scheme 
with CTEQ6 PDFs is used. 

Considering the scale dependence in the transition from LO to NLO at the LHC (see 
FigurefTTI). only a modest reduction is observed if W-boson pairs in association with two hard 
jets are taken into account. This large residual scale dependence is mainly due to the qg chan- 
nels, followed by those contributions with two valence quarks in the initial state that are present 
in the real corrections, but not at LO. The scale dependence can be significantly suppressed 
upon applying the veto of having "no 2 nd separable jet". The relevance of a jet veto in order 
to suppress the scale dependence at NLO was also realized for genuine W-pair production at 
hadron colliders [62]. A reduction of the difference between the two curves, which represents — 
as mentioned in the beginning of this section — the contribution of genuine WW+2jets events, 
is also achieved by increasing the value of the cut on pxjeb which is illustrated in Figure[T2~1for 
WW+jet production with pT,jet,cut > lOOGeV in the LHC setup and pT,jet,cut > 50GeV for the 
Tevatron. Explicit numbers for the cross sections — both with and without leptonic W decays — 
are collected in Tables [Q-IH 

In general, the influence of the restriction on genuine WW+jet production via the described 
jet veto is not that large in the Tevatron setup (see Figure (TT]). This can be understood from 
the lower CM energy at the Tevatron: The energy for producing a second hard jet is available 
less frequently here, so that a stronger suppression of WW+2jets events is obtained. The fact 
that the difference between the two NLO observables strongly decreases when going to higher 
/?T,jet,cut values — as shown on the right-hand side of Figure [12] — confirms this interpretation. 

In Figure [13] the respective plots are shown with the W decays included via the iNWA and 
with further cuts applied according to Eqs. (15.81) - (15.101) . Qualitatively, the features discussed 
above can be observed as well. 
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pp -> W + W- +jet + X @ 14TeV 


M — A*ren — A'fact 


o"Lo[pb] 


0"NLO,excl[pb] 


C"NLO,incl[pb] 


0.5 M W 
1M W 

2M W 


27.17638(89) 
23.97398(79) 
21.26027(70) 


22.216(15) 
22.606(12) 
22.5646(97) 


34.641(14) 
31.970(11) 
29.7959(89) 


pp -> W + (^ v e e + )W-(^ fi-v^+jet + X @ 14TeV 


M — A*ren — A'fact 


o-Lotfb] 


0"NLO,excl[ft>] 


0"NLO,incl[fb] 


0.5 M W 

1M W 

2M W 


118.849(11) 
104.9482(94) 
93.1789(83) 


89.227(99) 
93.789(76) 
95.493(60) 


148.826(68) 
138.808(56) 
130.401(57) 



Table 1 : WW+jet cross sections at the LHC, calculated in the five-flavour scheme with CTEQ6 
PDFs. The cut on the transverse momentum of the jet is set to 50GeV. In the lower part of the 
table where the decays of the W bosons are included further cuts according to Eqs. (15.81) - (15.101) 
are applied. 



pp -> W + W^ +jet + X @ 14TeV 


M — f^ren — A'fact 


o-Lotpb] 


0"NLO,excl[pb] 


0"NLO,incl[pb] 


0.5 M W 
1M W 

2M W 


12.01935(44) 
10.37216(39) 
9.02371(34) 


11.8943(76) 
11.4827(57) 
11.0114(46) 


16.2377(73) 
14.6889(57) 
13.4405(44) 


pp -> W + (^ v e e + )W"(^ fTvJ+iet+X @ 14TeV 


M — A'ren — /^fact 


o-Lotfb] 


C"NLO,excl[fb] 


0"NLO,incl[fb] 


0.5 M W 
1M W 

2M W 


61.6020(66) 
53.2170(57) 
46.3585(50) 


57.471(40) 
56.734(47) 
55.192(38) 


79.722(45) 
73.199(39) 
67.680(28) 



Table 2: As in TableCQ but with p T ,jet,cut = 100GeV. 
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pp -> W + W" + jet + X @ 1 .96 TeV 


M - Mien - /^fact 


& LOLP D J 


C^NLO.exclLPDJ 


& NLO,inclLP D J 


0.5 M W 
1M W 

2M W 


2.88967(29) 
2.30259(24) 
1.87503(20) 


2.6343(32) 
2.5761(24) 
2.4208(18) 


3.1417(30) 
2.9257(21) 
2.6688(16) 


pp^ W + (^ v e e + )W-(^yu-v^)+jet + X @ 1.96 TeV 


M - Mien - Mfact 


OLotfb] 


0"NLO,excl[fb] 


0"NLO,incl[fb] 


0.5 M W 

1M W 
2M W 


13.7407(21) 
10.9245(17) 
8.8789(14) 


12.2366(88) 
12.0582(72) 
11.3695(68) 


14.8002(91) 
13.8052(73) 
12.6152(68) 



Table 3: As in Tabled! but at the Tevatron with /?T,jet,cut = 20GeV. 



pp -> W+W" + jet + X @ 1 .96 TeV 


M — Mien — Mfact 


£TL0[fb] 


0"NLO,excl[fb] 


0"NLO,incl[fb] 


0.5 M W 

1M W 

2M W 


0.843556(80) 
0.657787(64) 
0.525619(52) 


0.76867(81) 
0.74241(58) 
0.68820(44) 


0.83508(80) 
0.78664(55) 
0.71884(44) 


pp^ W + (^ v e e + )W-(^ m~v^)+ jet + X @ 1.96 TeV 


M - Mien - Mfact 


cr LO [fb] 


0"NLO,excl[fb] 


0"NLO,incl[ft>] 


0.5 M W 

1M W 
2M W 


4.55105(67) 
3.54156(53) 
2.82493(43) 


4.0799(36) 
3.9594(27) 
3.6774(21) 


4.4484(36) 
4.2037(27) 
3.8466(21) 



Table 4: As in Table [H but at the Tevatron with /?T,jet,cut = 

50GeV. 



( ii) Dependence on the transverse-momentum cut on the jet 

To close the discussion of the NLO cross sections for WW+jet that are not sensitive to the 
decays of the W bosons, the dependence on the cut applied to the transverse momentum of 
the jet is considered here. In Figure [H the /?T,jet,cut dependence is shown for the LHC and 
the Tevatron setup, respectively. In order to introduce a measure for the scale uncertainties 
of the cross sections, bands are depicted that correspond to a variation of \x = ju ren = /if act by a 
factor 2 around the central scale. A crossing of the curves for differing scale values leads to 
vanishing band widths at some points, which is an artifact of how the results are depicted and 
should not be misinterpreted as a vanishing scale uncertainty. In the plots of Figure [141 the band 
corresponding to the more inclusive cross section is partially covered by the more exclusive one. 
Both curves are, however, shown here in order to confirm the statement that the effect of genuine 
WW+2jets events decreases with an increasing value of PT,jet,cut- This becomes manifest in the 
overlap of the two bands especially for large cut values at the Tevatron, but is also evident in 
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Figure 12: Scale dependence of the WW+jet cross sections with /u = yUf act = /u T en, where 
PTjet > lOOGeV is applied for the LHC, and PTjet > 50GeV for the Tevatron. 

the plots for the LHC setup from the convergence of the two NLO bands when going to larger 
values for px,jet,cut- Eventually, the PT,jet,cut variation reflects the behaviour discussed in the 
previous paragraph: For Tevatron, a considerable reduction of the scale uncertainty is achieved 
when going from LO to NLO, whereas this reduction is only mild for LHC unless WW+2jets 
events are vetoed. 

5.3 Results on differential cross sections 

Now we mainly focus on differential distributions in observables defined from the W decay 
products, i.e. the cuts on the decay leptons specified in Eqs. (15.8I) - (15.10I) are applied. 

5.3. 1 LO analysis of the different decay descriptions 

The integrated WW+jet cross section as defined above is not sensitive to the decays of the 
W bosons. An inclusion of leptonic decays into the calculation of WW+jet production can be 
performed in different ways. In this section, a comparison of LO results for the three strategies 
discussed in Section |4]is performed, which are a full amplitude calculation within the complex- 
mass scheme, the naive NWA, and its improved version iNWA that treats the W bosons as 
on-shell particles, but keeps spin correlations. The aim of this discussion is to find an adequate 
approximation in order to avoid performing the full amplitude calculation but still to obtain an 
appropriate description of the decays. 

Considering only integrated cross sections, both approximations reproduce the full results 
in the LHC setup at the expected accuracy of 0(Tv/My). More precisely, we find a deviation 
of about 3% for the naive NWA and of less than 1% for the improved version with /?T,jet,cut = 
50, 100 GeV. In the Tevatron setup, the naive NWA predicts cross sections deviating by roughly 
15% (9%) for p T jet,cut = 20 GeV (50 GeV), while the iNWA is good within 1%. The differences 
between the two approximations are due to the ignored spin correlations in the naive NWA, 
which lead to different distributions of the decay leptons over the phase space. These changes 
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Figure 13: Scale dependence of the WW+jet cross sections with W decays included and further 
cuts applied according to Eqs. (15.81) - (15.101) . The five-flavour scheme with CTEQ6 PDFs is 
used. In the LHC setup, the results are given for /?T,jet,cut = 50GeV (upper left plot) and for 
Px,jet,cut = lOOGeV (lower left plot). For the Tevatron we show results for /?x,jet,cut = 20GeV 
(upper right plot) and for px,jet,cut = 50GeV (lower right plot). 
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Figure 14: Variation of the WW+jet cross sections with pT,jet,cut at the LHC (l.h.s.) and at the 
Tevatron (r.h.s.): The bands correspond to a scale variation, again with jj, = fif act = /i ren , by a 
factor 2 around the central scale. The red and the green curves correspond to the different defi- 
nitions of the NLO observable, where the red one depicts the more inclusive quantity. The upper 
plots correspond to stable W's, the lower plots include leptonic W decays with the additional 
cuts of Eqs. d578b - (15.101) applied. 
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Figure 15: Comparison of different W-decay descriptions in the pj distributions of each of the 
two decay leptons. The LO cross sections are evaluated at p = fif act = p ren = M\y for the full 
amplitude calculation, the naive NWA, and the improved NWA. 



become manifest in the integrated cross sections if cuts on the decay products are applied. In 
general, both the naive and the improved version of the NWA give the same integrated cross 
sections if no leptonic cuts are applied. This is due to the fact that taking the spin correlations 
into account only causes a (W-polarization-dependent) redistribution of the momenta of the 
decay leptons, which does not influence the result as long as the full decay phase spaces are 
integrated over. 

Distributions in specific variables are, however, strongly affected, which can be seen from 
the differential cross sections for the LHC given in Figures [151 EH and [T7J Only the dis- 
tributions for decay leptons are shown there, because these are, naturally, mainly affected by 
the different decay descriptions. Considering the distributions of transverse momentum pi and 
pseudo-rapidity j] of each of the two decay leptons, the improved NWA delivers a very accurate 
reproduction of the full calculation, whereas the naive version deviates by up to 15% in some 
phase-space regions. The distributions of the angles between the two leptons — <p denotes the an- 
gle in the transverse plane and cos 6 the cosine of the angle between the two leptons — resulting 
from the full amplitude calculation are also in general reproduced more precisely by the im- 
proved NWA. At the Tevatron the situation is quite similar, which can be read off Figures [T8] 
and [T9j While the improved NWA reproduces the full result quite well, the predictions of the 
naive NWA deviate quite strongly — in particular in case of the pseudo-rapidity distributions of 
the leptons. 

The analysis of this section justifies the application of the improved NWA for the NLO QCD 
calculations to pp/pp — » WW + jet + X with leptonic decays. 
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Figure 16: As in Figure[T5l but for the rj distributions of the two decay leptons. 
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Figure 17: As in Figure [T51 but for the cosine of the opening angle, cos#, and the angle <p in the 
transverse plane between the two leptons. 
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Figure 18: Comparison of different W-decay descriptions in the pi and rj distributions of each 
of the two decay leptons. The LO cross sections are evaluated at fi = /if act = /i ren = Mw for the 
full amplitude calculation, the naive NWA, and the improved NWA at the Tevatron. 




Figure 19: As in Figure [HI but for the cosine of the opening angle, cos#, and the angle ip in the 
transverse plane between the two leptons. 
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Figure 20: Differential cross sections for WW+jet with decays included in the improved NWA 
at the LHC: The LO and NLO distributions are shown for /i = //f act = // ren = Mw. The distri- 
butions for the transverse momenta pj of the jet and of the decay leptons, and for the missing 
transverse momentum /?T,miss are depicted. The bands in the ^-factors are precisely defined in 
the text. 
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5.3.2 Differential NLO cross sections at the LHC 

A survey of distributions for differential WW+jet cross sections with leptonic W decays is 
given in Figures [20]— [22] for the LHC setup with px,jet,cut = 50GeV and p. = p ren = yUfact = Mw- 
The relative size of the NLO corrections is represented in terms of a ^-factor. The depicted 
bands correspond to a variation of the scale p by a factor of 2 in the NLO quantities only, i.e. 
we show <t NL o(//)/ctloC^%) with 0.5Mw <p< 2M\y. 

The jet distributions are understood as distributions of the hardest jet. Figure [201 provides 
transverse-momentum distributions of the hadronic jet, of the missing transverse momentum 
due to (anti-)neutrinos leaving the detector undetected (upper plots), and of the two charged 
decay leptons (lower plots). For all px distributions, a tendency of the more exclusive NLO 
cross section to decrease faster than the LO cross section when going to higher values is evident. 
To a large extent this is due to the fact that a fixed value is used for the renormalization scale. 
Instead, the transverse momentum of the jet seems a more appropriate scale in the high-px tail, 
since the only arising strong coupling concerns this jet. Using a s (Mw) in the LO calculation 
overestimates the contributions for large Px.jet values due to the ignored decrease of the QCD 
coupling. In other words, the large Pxjet introduces an additional scale which is responsible 
for a large logarithmic enhancement. These corrections may be absorbed into the running of 
the strong coupling constant (and thus resummed to all orders) by using a more appropriate 
renormalization scale, i.e. a scale that is set through pxjet- Similar arguments can explain part 
of the behaviour at small pj observed in the respective plots of Figure [201 but of course the 
theoretical description in regions of small transverse momenta should eventually be improved 
by dedicated resummations of higher-order QCD corrections. 

Since the leptonic transverse momenta are connected to Px.jet by momentum conserva- 
tion, the same argumentation also holds for these distributions. The effects are, however, 
much stronger for the jet compared to the leptons due to the tendency of coloured particles to 
collinearly radiate further QCD partons. For the more inclusive NLO observable, this effect is 
overcompensated by genuine WW+2jets events which are — being actually LO contributions — 
also influenced by this effect. In the pj distributions for the two decay leptons, a slight ten- 
dency to higher pj can be found for the lepton. That a difference between lepton and antilepton 
arises at all can be understood from the defined order of the W bosons coupling to the fermion 
chain, because /-channel-like emissions of the W bosons cause large contributions. Whereas 
a £-channel-like emission of a W + can only arise in case of an incoming up quark or down 
antiquark, an analogous emission of a W" always stems from an incoming down quark or up 
antiquark. The antilepton always results from the W + decay and the lepton from the W" decay. 
Due to the difference of the involved PDFs, the antilepton prefers lower transverse momenta 
compared to the lepton. 

An important impact of this effect can also be seen in the upper plots of Figure I2T1 where 
distributions of the pseudo-rapidities of the leptons are depicted. Here, the antilepton shows a 
slightly larger tendency to small angles against the beam axes than the lepton. Naturally, due to 
the symmetric hadronic initial state in proton-proton collisions the pseudo-rapidity distributions 
are symmetric. The size of the relative NLO corrections turns out to be nearly independent of 
these quantities — well described by a constant .K-factor. The symmetry property still holds for 
the pseudo-rapidity of the hadronic jet, which is shown in the lower plot of Figure [2T] The NLO 
corrections to this quantity, however, increase for larger pseudo-rapidities. These large positive 
corrections in the region of large pseudo-rapidities can be understood by a simple statistical 



29 



effect: NLO corrections in general redistribute events via radiation. Since most events show 
jets in the region of lower pseudo-rapidities, it is just more likely to redistribute events from the 
low-pseudo-rapidity region to the high-pseudo-rapidity region than vice versa. Apart from this 
small phase-space region, which is negligible anyway if the total amount of events is considered, 
the NLO corrections turn out to be nearly independent of the jet pseudo-rapidity as well. 

The invariant-mass distribution of the charged leptons, which is shown in the lower right 
plot of Figure [2H depicts the typical behaviour of two particles that are not resonantly produced: 
while the steep increase at low invarint masses is mainly due to the lep ton- separation cut and 
the transverse-momentum cuts on the individual particles, the decrease to higher values reflects 
the behaviour of all variables involving a new energy scale, according to the dependence of 
hadronic cross sections on the partonic centre-of-mass energy. 

The two distributions in Figure [22] depict the angular correlations between the two decay 
leptons. The angle between the two leptons in the transverse plane is represented by <p, which 
is a quantity invariant under boosts along the beam direction. By 6 the opening angle between 
these two leptons in the laboratory frame is denoted. Considering the <p distribution, i.e. ignor- 
ing the boost effect along the beam axes, the two charged leptons turn out to fly preferentially 
into opposite directions. This is not surprising, since momentum conservation forces the two 
W bosons at least to show a tendency to opposite directions which is mediated to their decay 
products by boost effects. The angle between the leptons is, however, important for the dis- 
tinction of the background process WW+jet from the signal process H(— > WW*)+jet in Higgs 
searches. This is due to the fact that the decay leptons of a W-boson pair arising from the decay 
of a scalar Higgs particle show the — on the first view non-intuitive — tendency to fly into the 
same direction. This property results from the spin correlation of the W + W~ system, which is 
discussed in detail in Ref. [63]: Since, in the rest frame of the Higgs boson, the helicities of 
the W bosons are correlated, and only a left-handed charged lepton and a right-handed charged 
antilepton can arise from the W decays, their emission in the same direction is favoured. This 
correlation effect is, of course, smeared in H(— > WW*)+jet, since the Higgs boson is in general 
boosted, but a remainder of the effect should still be measurable. 

For WW+jet production, however, the helicities of the two W bosons are not correlated in 
this way, as depicted in Figure [221 but tend to opposite directions. In the cos 6? distribution, this 
preference is overcompensated by the boost effect along the beam axes, leading to a tendency 
in direction of small opening angles. The dependence of the size of NLO corrections on both 
angles turns out to be of the order of 10-20%. 

5.3.3 Differential NLO cross sections at the Tevatron 

Considering the same quantities for WW+jet production in the Tevatron setup, most of the 
effects can be explained by the fact that proton-antiproton collisions take place here instead 
of proton-proton collisions. Besides, the lower CM energy compared to the LHC plays an 
important role. We start again with the discussion of transverse-momentum distributions of 
the hard jet and the decay leptons, and of the missing transverse momentum, which are shown 
in Figure [23) Here, the behaviour of the more exclusive NLO observables can in principle be 
explained in the same way as for LHC: the LO results overestimate the cross section at high 
scales and underestimate it at low scales due to the fixed renormalization scale used in the 
calculation. As already observed when considering integrated cross sections, the difference 
between the two NLO observables is quite small at the Tevatron, which can be understood 
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Figure 21: As in Figure [201 but for the pseudo-rapidity rj of the charged decay leptons (upper 
plots), the pseudo-rapidity of the jet (lower left plot), and the invariant mass of the charged 
leptons (lower right plot). 
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from the smaller CM energy: In most cases, not enough energy is available for the production 
of a second hard jet. Therefore, the LO contributions of pp — » WW + 2jets+X contained in 
the more inclusive NLO observable only weaken the effect of negative corrections at large pj, 
but do not overcompensate it as for the LHC setup. The pj distributions of the two decay 
leptons are identical up to numerical fluctuations. This results from the fact that the hadronic 
process pp — » WW+jet + X and the underlying model (SM without CP- violating phases) are 
invariant under CP transformation. Having not introduced any such CP violating contributions, 
we thus expect that the final state is also even under CP transformation. Thus, in contrast to 
proton-proton collisions, PDF effects do not cause differences between the pi distributions of 
the leptons, since the quark PDFs of the proton equal the antiquark PDFs of the antiproton. 

The distributions of the pseudo-rapidities j] of the leptons, which are depicted in Figure l24l 
are not symmetric with respect to J] = due to the asymmetric hadronic initial state. Instead, the 
distributions of the lepton and the antilepton are identical if one of the distributions is mirrored 
around the r\ = axes. In the Tevatron setup, the positive beam axes corresponds to the direction 
of the proton beam, and the negative axes to the antiproton beam. Correspondingly, positive 77 
values describe momenta tending in the direction of the proton beam and vice versa. 

Considering the non-symmetric pseudo-rapidity distributions, a tendency of the antilepton 
to the positive beam direction and, correspondingly, of the lepton to the negative beam direction 
is evident. To explain this, the argument of /-channel-like emission of W bosons can again be 
applied. As discussed in the LHC case, the W + boson, and consequently the antilepton, can only 
be emitted in this way from an up quark or a down antiquark. Since the qq channels dominate 
at the Tevatron, only these have to be taken into account for a qualitative discussion. On the 
average, these subprocesses do not produce boosted events in a distinguished direction due to 
the fact that the quark PDFs in the proton are equal to the respective antiquark PDFs in the 
antiproton. The contributions u p (pi)Up(/?2) and d p (pi)dp(p2) dominate over the contributions 
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Figure 23: Differential cross sections for WW+jet with decays included in the improved NWA 
at the Tevatron: The LO and NLO distributions are shown for fi = /if act = /i ren = M\y. The dis- 
tributions for the transverse momenta pj of the jet and of the decay leptons, and for the missing 
transverse momentum /?T,miss are depicted. The bands in the ^-factors are defined as in the 
LHC case. 
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Figure 24: As in Figure [231 but for the pseudo-rapidity rj of the charged decay leptons (upper 
plots), the pseudo-rapidity of the jet (lower left plot), and the invariant mass of the charged 
leptons (lower right plot). 
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Figure 25: As in Figure [2U but only the LO distributions for the pseudo-rapidities J] of the 
lepton and the antilepton are depicted. In addition, the two dominating partonic contributions 
uu and dd to the hadronic cross section are shown. 

u p (pi)Up(p2) and d p (pi)dp(/?2)- Due to the effect of /-channel-like W + emission, the first tends 
to antilepton emission in the positive beam direction, the latter to the negative direction, as 
depicted in Figure [251 Since the uu channel exceeds the dd channel roughly by a factor 4 in 
total, a tendency to the positive beam direction results. The same argumentation holds, mutatis 
mutandis, for the lepton and the W" boson. 

The pseudo-rapidity of the jet, which is depicted in the lower left plot of Figure [231 is sym- 
metric with respect to rj = 0. This is due to the fact that no distinction can be made between 
hadronic jets arising from gluons, quarks, or antiquarks. Therefore, the sum over all contri- 
butions, which are CP symmetric in pairs, yields a symmetric distribution. As in the case of 
LHC, the dependence of the NLO corrections on the pseudo-rapidities of the jet and the lep- 
tons is moderate. Considering the absolute pseudo-rapidity dependence of the cross sections, a 
tendency to events that are not strongly boosted in the direction of the beam axes is observed. 
This can be understood from the fact that the dominating partonic channels with initial states of 
valence quark and valence antiquark are not strongly boosted in general. 

The behaviour of the invariant-mass distribution of the two charged leptons, which is de- 
picted in the lower right plot of Figure [24l can be explained analogously to the LHC case. 

Finally, the angle correlations between the two decay leptons are considered. The angle 
in the transverse plane is again labelled by ip, and 9 is the opening angle between the two 
leptons. The corresponding distributions are depicted in Figure[26l These angle correlations are 
very important also at Tevatron to distinguish the background process WW+jet from the signal 
process H(— » WW*)+jet. Details on this topic are given in the discussion of the respective LHC 
distributions. Also in the Tevatron setup, the two leptons tend to fly into opposite directions in 
the transverse plane. In contrast to the situation at the LHC, this tendency is still observed in the 
cos 9 distribution of the opening angle 9, i.e. the effect is not overcompensated by boost effects. 
This is again understood from the fact that no tendency to strongly boosted events arises in the 
dominant partonic channels at the Tevatron. 
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Figure 26: As in Figure [231 but for the angle <p in the transverse plane and the cosine of the 
opening angle between the two leptons. 



5.4 Benchmark results and comparison to results of other groups 

Independently of our calculation described in Ref. [9], two further groups have performed 
NLO QCD calculations on WW+jet production. The results of one of these groups are given 
in Ref. [10], while those of the second group are not published yet [11]. In Ref. [3], a tuned 
comparison of the three calculations is provided for the integrated LO cross section, which agree 
again within statistical errors, and of the virtual corrections at one specific phase-space point. 
There, the renormalized matrix elements are subdivided into bosonic and fermionic corrections 
with only the two light generations included in closed quark loops. The results are given in terms 
of the coefficients to the poles jj and i, and the finite part, i.e. in a way that is independent of 
the specific method for cancelling the infrared singularities @ A comparison of results for the 
third-generation loops has not been performed so far. Agreement is achieved between the results 
of the three groups at an accuracy level comparable to our internal checks. 

6 Conclusions 

The production of W-boson pairs in association with a hard jet is an important source for 
background to Higgs and new-physics searches both at the Tevatron and at the LHC. A proper 

4 Unfortunately an important detail was not explicitly stated in Ref. [10]: In the dipole subtraction method the 
IR divergences are usually factorized from D-dimenional LO structures ("conventional dimensional regulariza- 
tion"). In the comparison shown in Ref. [10], however, absolute numbers on the corrections to the squared matrix 
elements are given where the D-dimenional LO structures are evaluated in four dimensions (without expanding 
around D = 4). This procedure is legitimate, because the same replacement of D-dimensional by four-dimensional 
LO structures is done in the IR-divergent part of the real corrections. The numbers given in Ref. [10] are, thus, 
equivalent to results in the '"t Hooft-Veltman scheme" for dimensional regularization where only particles in loops 
or in singular splittings are extended to D dimensions. 
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theoretical prediction for this process requires at least the inclusion of perturbative QCD cor- 
rections at the next-to-leading order. 

Continuing previous work, we have described in detail the calculation of these corrections, 
which takes into account leptonic decays of the W bosons. To this end, we employ an improved 
narrow-width approximation that treats the W bosons as on-shell particles, but keeps the W- 
spin information. While the naive narrow-width approximation, which neglects the information 
on the W spins, is only good within 5-10% in differential distributions, the improved version 
roughly reaches percent accuracy. This is the result of a comparison of leading-order predictions 
in these two approximations with results of a calculation fully based on off-shell W bosons. 

The treatment of bottom quarks in the initial or final states deserves particular attention. 
Numerically such contributions only play a significant role if the top-quark propagator present 
in this case becomes resonant. The contributions are thus essentially the off-shell continua- 
tions of W + t, W"t, or tt production including the subsequent decay of the top-quarks. These 
subprocesses, however, should not be counted as part of WW+jet production. Therefore, we 
have excluded contributions from external bottom quarks. The reliability of our procedure has 
been verified in a comparison of next-to-leading-order results obtained in two different schemes, 
where one is based on four, another on five active quark flavours in the proton. 

Our detailed discussion of numerical results shows that the QCD corrections stabilize the 
leading-order prediction for the WW+jet cross section considerably with respect to a variation 
of the factorization and renormalization scales which we identify with each other. At the LHC, 
this stabilization of the prediction, however, requires a veto on a second hard jet. Otherwise 
the production of final states with WW+2jets, which yields a leading-order component of the 
next-to-leading-order correction, introduces again a large scale dependence. As far as the dif- 
ferential distributions are concerned the corrections are typically of the order of about 25%. For 
a remarkable number of distributions the ^-factor is only mildly dependent on the kinematical 
region. At the LHC the 77 distributions in the dominant region and also the distributions in the 
angles between the two charged leptons have an almost constant ^-factor of about 1.3 (inclu- 
sive cross-section definition); for the exclusive cross-section definition the corrections are even 
smaller and rather close to 1. The pi spectra, on the other hand, show a much more phase- 
space-dependent ^-factor with the exclusive cross-section definition showing an even larger 
dependence than the inclusive one. This is not surprising since the pj introduces an additional 
scale which could introduce potentially large logarithms which are badly treated by a constant 
renormalization scale. At the Tevatron our findings are similar. Again the rj and angular dis- 
tributions receive corrections in form of an almost constant ^-factor of about 1.3 (exclusive 
definition). The corrections for the exclusive cross-section definition are again smaller than for 
the inclusive definition. In case of the pj spectra we observe again a phase-space-dependent 
i^-factor. We note that the almost constant A"-factor which holds for a remarkable number of 
distributions has also been observed by Campbell, Ellis and Zanderighi [10]. 

The QCD corrections to the related processes of ZZ+jet and WZ+jet production can be 
obtained in an analogous way as presented here for WW+jet. The corresponding calculations 
are in progress. 
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Figure 27: Schematic diagram for the process ab — > VX — » ff'X. 
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Appendix 

A Decaying on-shell gauge bosons 



In this appendix we describe a simple method to include gauge-boson decays with the cor- 
rect spin correlations in processes for which the helicity amplitudes for on-shell gauge bosons 
are already known. Specifically we consider the situation illustrated in Figure [27] where the 
vector boson V with momentum k decays into the pair f(k\)f'(k.2) of massless fermions with 
momenta k\p_. We start with the expression for the cross section 



1 



ab-+VX^>ff'X 



2s a b 



dO 



ffx 



M 



ab->VX->ff'X 



(A.l) 



for the full reaction ab — » VX — » ff'X integrated over its phase space Oyy> x at tne CM energy 
■\fsab~. Of course, not only diagrams with a V resonance, as shown in Figure [27] contribute to 
the full amplitude M^^x^ffx- We are, however, interested in the resonant diagrams which 
factorize according to 



M 



(res) 

ab^VX^ff'X 



Tu(k) 



■1 



k 2 -M\ + iM v r v 



f(h,k 2 ), 



(A.2) 



where T^(k) describes the production of an off-shell V and j^{k\,k2) is the current of the V 
decay. For small decay widths Ty of V the phase-space integral (IA.1I) is dominated by the 
resonant diagrams with momenta near the mass shell of V. In order to extract the resonant 
terms, we factorize the phase space into V production and decay as follows, 



dO 



1 



ffx 



\k 2 -Ml + iM v r v | 2 
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dk 2 1 



In \k 2 -M 2 +iM v r v | 2 



d(S>Yx(k) I d® f f,(ki,k 2 ) 



1 



2M v r v 



dOvz(fe) / dO^i.fe)--, (A.3) 



where the hats over the momenta indicate that the on-shell condition k = (k\ + k 2 ) = My is ful- 
filled. On resonance we can decompose the contraction T^j^ into helicity amplitudes Mab^vx 
and Af v _»//' for the V production and decay upon inserting the completeness relation for the V 
polarization vectors s^(A), 

A=0,±l M V J 

= J2 M ab ^ YX (X)M Y ^ ff -,(A) (A.4) 

A=0,±l 

where we have used current conservation k v j v (k\,k2) = (for massless fermions) and identified 
Mab^vxW = T M (k)4(AT, M y ^ ff ,(A) = E Y (A)j v (k u k 2 ). (A.5) 
In summary, the cross section in resonance approximation reads 

O-ab^VX^ff'X rvTo Y~ I d °Vx(^) Yl M ab^Vx{X)*Mab^Vx(X) 

Sab ■' X,M=%±\ 

X 1m\^ I MtfitutoMy^AyMv^W. (A.6) 

If the phase space of the decay fermions is integrated over completely, rotational invariance 
implies that 

^- J d® ff 4iMM Y ^ ff ,(AyM Y ^ ff -,(A) = T y ^ ff -,6 AU , (A.7) 

where Ty^^y, is the partial decay width. Inserting this result into (IA.6I ), yields the usual narrow- 
width approximation for the cross section, 

o-ab^Yx^ffx f~ — I dQ V x(&) Yl \M ab ^vxW\ 2 BR y ^ ff , 

LSah J A=0,±l 

= o- ab ^vxBRy^ f f,, (A.8) 

with (Tai^Yx denoting the production cross section for unpolarized vector bosons V and BR v ^yy, 
r v ^yy,/rV denoting the branching ratio for the decay of V into the //' pair. 

If the decay fermions are not integrated over the full phase space ®ff>, i.e. if cuts on the 
fermions are applied or if distributions in the fermion kinematics are considered, it is convenient 
to introduce the "decay correlation matrix" 

A Y ^ ff '(h,k2) = 2MvIV AW/'^AW/'W. (A.9) 
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so that the full cross section takes the form 



ah- 



Sab J X,A'=0,±l 



xBR 



y^ff 



dO ff -,0q,fc 2 ) y _> f f, . . 



(A.10) 



where Oyy, is the volume of the //' phase space J dOff,(ki,k2). The matrix A v ^^ is widely 

independent of the production process. Only the V helicity states entering A v ^ f must be the 
same as in the production. The explicit calculation of A v_> ^ is conveniently performed in the 
rest frame of V. Defining the momentum and the polarization vectors of V in the CM frame 2 
of ab b>H 



P = Ey{l,pyCOS(pysm0y,Pysm(pysm6y,PyCOs6y), fiv = \/l-My/£y, 



£^,(±1) = — ^(O,-cos0vcos^v±isin^v>-sin0vcos^v + icos^v ? sin0v) ? 
V2 



£y(0) = 7v(y0v,cos0vsin0v,sin0vsin0v,cos0v), 7v = 1/ \/l -/?v = Ey/My 



(A.ll) 



we transform them into the rest frame Z of V with the Lorentz transformation matrix 



(A"v) 



/ yv -jSvrv\ /l o \ 

1 cos0vcos#v sin^ycos^v -sin0v 

1 -sin0v cos^v 

\-fivyy yv / V0 cos0vsin#v sin0vsin#v cos#v / 

(A.12) 

which is factorized into a rotation in S (rotating k to the z axis) and a boost along the new V 
direction. The explicit results for the considered vectors IT = A^ v a y in E are 



fr = (M y , 0,0,0), e£(±l) 



j+i0V 



■(o,-i, +i,o), 4(0) = (o,o,o,d 



(A.13) 



In 2 the decay momenta ki^ are distributed isotropically and parametrized by 

% II ~ ~ ~ ~ ~ A / i X "X J. 

K[ = — (I,cos0isin0i,sin0isin0i,cos0i), K 1 =kr-k\. (A. 14) 

Equipped with these kinematical definitions we evaluate the decay matrix elements 

My^ ff ,(A) = egy'fftUfihMA^a-Vfih), (A. 15) 

where we made the chirality cr = + explicit (which coincides with the sign of the helicity of 
f) and used an obvious notation for the Dirac spinors Uf, vp and the chirality projectors a>± = 



5 For e^j(A) this choice follows the phase convention of Ref. [43]. 
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(1 +75)/2. The chiral couplings gy//' are defined as in (14.51) . Ordering the A values according 
to (+,0, -), the decay amplitudes can be written as 

M^ ff ,(A) = eg$ ff ,Mv (+^e^'-^)(l + cos§ 1 ),sin0 1 ,-^e- 1 *-^)(i - COS 0i)) , 
M-^ ffl (X) = egy ff M Y (-^= e ^'-^)(l - cos X ), sin 0i, +^=6-^-^(1 + cos#i)) . 

(A.16) 

Inserting these amplitudes into (IA.9I ) and summing over fermion helicities, we obtain for the 
(hermitian) decay correlation matrix 



A v -//' = 3c + yff , 



( ±(1 +COS6>i) 2 * * \ 

cos £i) sin ±sin 2 #i * 

2i(^!-0 v ) . 2n z<fa-<h) n \ ■ n 1 /1 n \2 I 

\ - £ — 4 — sin z ^i ~ 9a /9 (l-cos6>i)sm6>i ±(l-cos0i)v 

1/1 mt ,2 \2 y 



f ^(l-COS0i) Z * * 

\ -^W^ ^(l+cos^sin^ id+cos^) 2 / 

(A. 17) 



where we used the shorthand 



4// ' = ^//') 2 + (^v//') 2 ' C v//' +c v//' = L (A. 18) 



From this result it is obvious that the full integration of the decay phase space yields 

V/' 

as it should be for an isotropic decay in the rest frame of V. 



J — ^ A v ^ (fc.fe) = y ^ A v ^ (fc.fe) = 1, (A.19) 



B Benchmark numbers for the virtual corrections 

In order to facilitate a comparison to our calculation, we provide explicit numbers on the 
squared LO amplitude and the corresponding virtual corrections for a single non-exceptional 
phase-space point. The set of momenta for ab — » W + (— > v e e + )W~(— > // _ v yU )c with the explicit 
partonic reactions qq — > W + W"g, qg — * W + W"q, and gq — > W + W"q is chosen as 

p a = (250,0,0,250), 

Pb = (250,0,0,-250), 

p Ve = (49.43764668100422,-18.32747340442861,-42.74042990560614,16.77618190917520), 

p e + = (175.4405115231198,-125.7805594116281,-26.45739027367693,119.4095073999478), 

Pyr = (32.78781503058663,-29.71117802999691,7.333312535152912,-11.76899493633845), 

Pv M = (81.55218598686336,45.19373677627630,61.86450764413016,-27.94759183511103), 

p c = (160.7818407784260,128.6254740697773,0,-96.46910253767348), (A.l) 
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with the obvious notation p = (p° , p l , p 2 , p 3 ) and all the components given in GeV. Note that 
the presented results are calculated using the improved NWA, i.e. with the W bosons on-shell, 
but with their spin information kept. 

Factoring out the couplings for the Born amplitude, we define 

\jf E \^o\ 2 = e A gla Q , (A.2) 

C spin, colour 

where 

IMloI 2 = M ah ^ w+w - c (A' + ,A'_TM ab ^ w+w - c (A + ,A-) (A3) 

A+,A' + ,A-,A'_=0,±1 

with the definitions of App.|Al The factor 1 /N c is due to the average over the incoming colour. 
For the channel qq, qg, gq we have N c = 9,24,24, respectively. 

The finite remainder of the virtual amplitudes after renormalization and addition of the I 
operator defines the relative correction cq according to 

-— 2Re ( (Mvirt + M ct + Mj-op) *M L o) =e 4 g 2 s a c , (A.4) 

spin, colour 

where A4 v i rt , Ai c u an d A4j_ p are defined in analogy to (IA.3I ). 

The results of our independent calculations at the sample phase-space point (IA.1I) are col- 
lected in Table |5] in terms of ao and Co, namely for the six different partonic insertions for a, b, 
and c. 

For the LO amplitudes we find an agreement of at least 14 digits — as expected from a 
calculation using 64bit double precision. For the finite remainder of the one-loop corrections 
we find an agreement of at least 8 digits. 

C Tables for histograms 

In this appendix we present the tables corresponding to the differential distributions pre- 
sented in Section |5.3.21 For each distribution, we list all predictions for the LO cross section 
and for both definitions of NLO cross sections — namely NLO (exel) and NLO (incl) — with the 
scale choice /j. = /i ren = /"fact = Mw- The given errors result from the Monte Carlo integration. 
The bin is specified by its central value. The bin width — which we chose constant for the entire 
histogram — is obtained from the distance of two neighboring bin positions. 
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[GeV] 



dcr L0 
d/'r 



Pb 
GeV 



,excl 
^PTjniss 



Pb 
GeV 



27.5 

32.5 
37.5 
42.5 
47.5 
52.5 
57.5 
62.5 
67.5 
72.5 
77.5 
82.5 
87.5 
92.5 
97.5 
102.5 
107.5 
112.5 
117.5 
122.5 
127.5 
132.5 
137.5 
142.5 
147.5 



1.10787(34) 
1.21978(36) 
1.28976(37) 
1.32096(37) 
1.31712(37) 
1.28158(36) 
1.22095(35) 
1.14291(33) 
1.06101(31) 
0.97356(30) 
0.88876(28) 
0.80514(26) 
0.72506(24) 
0.64875(23) 
0.57792(21) 
0.51221(20) 
0.45154(18) 
0.39710(17) 
0.34905(16) 
0.30907(15) 
0.27483(14) 
0.24595(13) 
0.22077(12) 
0.19854(11) 
0.17986(11) 



1.1683(35) 
1.2943(38) 
1.3484(33) 
1.3665(37) 
1.3687(38) 
1.3103(34) 
1.2395(30) 
1.1372(30) 
1.0611(28) 
0.9470(26) 
0.8481(24) 
0.7572(22) 
0.6524(30) 
0.5764(35) 
0.4942(25) 
0.4267(23) 
0.3666(19) 
0.3095(15) 
0.2663(18) 
0.2268(18) 
0.1925(17) 
0.1667(23) 
0.1474(24) 
0.1263(20) 
0.1137(14) 



do- nlo ,incl 
'//>! , 



Pb 

GeV 



1.5232(31) 
1.6798(31) 
1.7642(28) 
1.7964(31) 
1.8037(27) 
1.7481(27) 
1.6677(26) 
1.5551(26) 
1.4507(25) 
1.3188(23) 
1.1929(20) 
1.0743(18) 
0.9513(25) 
0.8452(23) 
0.7419(20) 
0.6576(20) 
0.5769(17) 
0.5058(13) 
0.4449(13) 
0.3939(13) 
0.3474(12) 
0.3096(19) 
0.2794(18) 
0.2487(12) 
0.2262(10) 





d<r LO 

dPTjet 


r P b " 

GeV 




do NLO.excl 
f 'PT,jet 


r P b " 

GeV 




do NLO.incl 
^PTjet 


r P b " 

GeV 




55 


1.47878(26) 


1.7903(23) 


1.8600(22) 


65 


1.19041(23) 


1.3622(21) 


1.5178(20) 


75 


0.98096(21) 


1.0595(20) 


1.2619(19) 


85 


0.82324(19) 


0.8511(20) 


1.0706(20) 


95 


0.70022(18) 


0.6875(19) 


0.9150(19) 


105 


0.60034(16) 


0.5648(18) 


0.7878(16) 


115 


0.51805(15) 


0.4642(21) 


0.6844(18) 


125 


0.45004(14) 


0.3871(27) 


0.5951(23) 


135 


0.39293(13) 


0.3342(25) 


0.5301(21) 


145 


0.34383(12) 


0.2711(14) 


0.4581(13) 


155 


0.30246(11) 


0.2306(17) 


0.4047(15) 


165 


0.26699(10) 


0.1952(22) 


0.3597(14) 


175 


0.236277(94) 


0.1721(20) 


0.3209(11) 


185 


0.209810(88) 


0.1448(10) 


0.28598(99) 


195 


0.186526(82) 


0.11976(96) 


0.25141(95) 


205 


0.166336(77) 


0.1070(10) 


0.23005(96) 


215 


0.148596(72) 


0.0883(12) 


0.20236(95) 


225 


0.133168(67) 


0.0796(13) 


0.1853(10) 


235 


0.119463(63) 


0.0666(11) 


0.1643(12) 


245 


0.107421(59) 


0.05777(98) 


0.1487(10) 


255 


0.096866(56) 


0.05028(86) 


0.13439(82) 


265 


0.087248(52) 


0.04629(82) 


0.12439(81) 


275 


0.078964(49) 


0.03675(67) 


0.11029(65) 


285 


0.071329(46) 


0.03121(58) 


0.09900(50) 


295 


0.064789(44) 


0.02906(56) 


0.09203(49) 





dcr LO 
dpTp- 


r P b " 

GeV 




do NLO.cxcl 


r P b " 

GeV 




do NLO.incl 

dPiti- 


r P b " 

GeV 






pt p+TGeVl 

y i ,e i * i 


do LO 
dPi fi * 


r P b " 

GeV 




d(T NLO.excl 
dPT,e+ 


r P b " 

GeV 




do NLO.incl 
dPT,e+ 


r P b l 

GeV 


27.5 


1.75438(39) 


1.7529(22) 


2.3634(20) 




27.5 


2.21854(59) 


2.1958(33) 


3.0506(27) 


32.5 


1.72719(38) 


1.7365(37) 


2.3444(34) 




32.5 


2.08923(52) 


2.0766(38) 


2.8827(28) 


37.5 


1.63647(37) 


1.6451(33) 


2.2317(31) 




37.5 


1.89828(46) 


1.8785(32) 


2.6195(30) 


42.5 


1.51599(36) 


1.5200(22) 


2.0711(20) 




42.5 


1.69226(41) 


1.6640(27) 


2.3308(25) 


47.5 


1.38745(35) 


1.3796(26) 


1.8914(21) 




47.5 


1.49414(37) 


1.4507(26) 


2.0407(23) 


52.5 


1.25864(34) 


1.2391(24) 


1.7103(22) 




52.5 


1.31341(34) 


1.2568(24) 


1.7847(22) 


57.5 


1.13443(32) 


1.0967(22) 


1.5288(20) 




57.5 


1.15140(31) 


1.0901(23) 


1.5565(22) 


62.5 


1.01762(31) 


0.9754(24) 


1.3693(20) 




62.5 


1.00791(28) 


0.9346(21) 


1.3457(21) 


67.5 


0.91166(29) 


0.8638(24) 


1.2242(19) 




67.5 


0.88205(25) 


0.8052(18) 


1.1693(18) 


72.5 


0.81508(28) 


0.7534(22) 


1.0821(19) 




72.5 


0.77219(23) 


0.6910(19) 


1.0169(17) 


77.5 


0.72755(26) 


0.6632(22) 


0.9638(19) 




77.5 


0.67635(21) 


0.5999(20) 


0.8876(16) 


82.5 


0.65001 (25) 


0.5885(21) 


0.8634(21) 




82.5 


0.59413(20) 


0.5165(20) 


0.7754(14) 


87.5 


0.58101(24) 


0.5157(18) 


0.7657(16) 




87.5 


0.52221(18) 


0.4458(16) 


0.6764(13) 


92.5 


0.51887(22) 


0.4477 (24) 


0.6805(19) 




92.5 


0.45992(17) 


0.3870(12) 


0.5926(11) 


97.5 


0.46488(21) 


0.3969(23) 


0.6081(18) 




97.5 


0.40599(15) 


0.3320(13) 


0.5189(12) 


102.5 


0.41717(20) 


0.3471(21) 


0.5405(18) 




102.5 


0.35941(14) 


0.2878(12) 


0.4557(10) 


107.5 


0.37503(19) 


0.3020(32) 


0.4837(26) 




107.5 


0.31924(13) 


0.2541(12) 


0.40510(99) 


112.5 


0.33748(18) 


0.2761 (29) 


0.4415(27) 




112.5 


0.28390(12) 


0.2169(11) 


0.35378(99) 


117.5 


0.30427(17) 


0.2373(29) 


0.3903(21) 




117.5 


0.25330(12) 


0.1941(11) 


0.31802(98) 


122.5 


0.27497(16) 


0.2180(31) 


0.3534(17) 




122.5 


0.22633(11) 


0.1665(10) 


0.27950(77) 


127.5 


0.24876(15) 


0.1870(16) 


0.3175(13) 




127.5 


0.20244(10) 


0.1481(11) 


0.24936(85) 


132.5 


0.22548(14) 


0.1671(13) 


0.2866(11) 




132.5 


0.181856(95) 


0.13194(98) 


0.22408(80) 


137.5 


0.20467(13) 


0.1502(13) 


0.2609(11) 




137.5 


0.163585(88) 


0.11121(72) 


0.19703(67) 


142.5 


0.18639(13) 


0.1313(13) 


0.2340(11) 




142.5 


0.147512(83) 


0.10190(72) 


0.17995(69) 


147.5 


0.16965(12) 


0.1176(13) 


0.2136(12) 




147.5 


0.133188(78) 


0.08883(82) 


0.16071(71) 



daro 
d V f r 



[pb] 



,excl 



[pb] 



^Q"NLO.incl 



[pb] 



-2.4 
-2.2 
-2 
-1.8 
-1.6 
-1.4 
-1.2 
-1 
-0.8 
-0.6 
-0.4 
-0.2 

0.2 
0.4 
0.6 
0.8 
1 

1.2 
1.4 
1.6 
1.8 
2 
2.2 
2.4 



11 

13 
15 
17 
19 
21 
22 
23 
24 
25 
26 
26 
26 
26 
26 
25 
24 
23 
22 
21 
19 
17 
15 
13 
11 



8284(58) 
7789(61) 
7736(63) 
7086(66) 
5103(68) 
1445(70) 
5806(71) 
8165(72) 
7995(72) 
5677(72) 
1199(72) 
.4622(73) 
.5492(72) 
.4622(73) 
.1199(72) 
5677(72) 
7995(72) 
8165 (72) 
5806(71) 
1445(70) 
5103(68) 
7086(66) 
7736(63) 
7789(61) 
8284(58) 



10.997(56) 
12.593(65) 
14.326(85) 
15.974(63) 
17.584(61) 
18.862(81) 
20.112(61) 
21.145(60) 
21.871(59) 
22.640(67) 
23.078(63) 
23.332(86) 
23.442(54) 
23.457(80) 
23.089(64) 
22.662(67) 
21.937(70) 
21.132(65) 
20.187(66) 
18.763(72) 
17.552(66) 
16.061(68) 
14.322(84) 
12.597(68) 
11.051(47) 



15.461(42) 
18.033(73) 
20.574(48) 
23.258(51) 
25.653(47) 
27.866(51) 
29.951(46) 
31.538(41) 
32.946(38) 
34.059(41) 
34.710(41) 
35.232(50) 
35.368(44) 
35.242(46) 
34.800(46) 
34.004(44) 
32.954(44) 
31.538(39) 
29.956(47) 
27.798(45) 
25.651(54) 
23.201 (64) 
20.606(45) 
18.044(44) 
15.473(37) 



77,,+ 


dr] e+ L P UJ 


ATNLO,exclr 

rf^ e + L P UJ 


dr] e+ L P UJ 


-2.4 


13.6439(77) 


12.462(66) 


17.890(69) 


-2.2 


15.4666(77) 


14.058(74) 


20.255(64) 


-2 


17.1745(76) 


15.362(55) 


22.547(42) 


-1.8 


18.7412(74) 


16.866(62) 


24.688(43) 


-1.6 


20.1160(71) 


17.900(69) 


26.546(45) 


-1.4 


21.3037(69) 


18.923(73) 


28.009(46) 


-1.2 


22.3357(67) 


19.83(11) 


29.498(51) 


-1 


23.1671(66) 


20.620(60) 


30.737(43) 


-0.8 


23.8414(64) 


21.168(77) 


31.588(46) 


-0.6 


24.3613(63) 


21.733(51) 


32.389(37) 


-0.4 


24.7380(62) 


22.134(64) 


32.918(37) 


-0.2 


24.9600(62) 


22.213(67) 


33.193(39) 





25.0380(62) 


22.397(68) 


33.450(62) 


0.2 


24.9600(62) 


22.253(68) 


33.187(55) 


0.4 


24.7380(62) 


22.121(65) 


32.998(43) 


0.6 


24.3613(63) 


21.697(53) 


32.355(42) 


0.8 


23.8414(64) 


21.189(65) 


31.653(43) 


1 


23.1671(66) 


20.550(56) 


30.653(42) 


1.2 


22.3357(67) 


19.716(90) 


29.496(56) 


1.4 


21.3037(69) 


18.996(73) 


28.069(49) 


1.6 


20.1160(71) 


17.926(64) 


26.477(50) 


1.8 


18.7412(74) 


16.845(66) 


24.663(50) 


2 


17.1745(76) 


15.363(57) 


22.498(41) 


2.2 


15.4666(77) 


13.978(67) 


20.250(51) 


2.4 


13.6439(77) 


12.439(67) 


17.922(47) 



m^+u- TGeVl 


do\n 
d>n c + t r 


r P b " 

GeV 




d<T NLO.cxcl 


r P b " 

GeV 




dT NLO.incl 

dm e+ll - 


r P b " 

GeV 




5 


0.013441(30) 


0.01412(22) 


0.01878(20) 


15 


0.139490(90) 


0.13426(44) 


0.18829(40) 


25 


0.24176(12) 


0.22708(55) 


0.32658(51) 


35 


0.32454(14) 


0.30439(67) 


0.43912(58) 


45 


0.40439(15) 


0.38108(76) 


0.54827(68) 


55 


0.50351(17) 


0.48153(84) 


0.68487(77) 


65 


0.60746(19) 


0.5910(10) 


0.82857(89) 


75 


0.68007(20) 


0.6630(11) 


0.92526(99) 


85 


0.70790(21) 


0.68676(98) 


0.95872(91) 


95 


0.69651(21) 


0.6692(12) 


0.9389(10) 


105 


0.65662(20) 


0.6262(12) 


0.8842(11) 


115 


0.60452(19) 


0.5627(11) 


0.8069(10) 


125 


0.54407(18) 


0.5029(11) 


0.72856(98) 


135 


0.48520(17) 


0.4442(11) 


0.64768(98) 


145 


0.42969(16) 


0.3878(10) 


0.57187(85) 


155 


0.37823(15) 


0.3324(10) 


0.49878(89) 


165 


0.33284(14) 


0.29197(93) 


0.44035(85) 


175 


0.29239(13) 


0.25088(80) 


0.38384(66) 


185 


0.25692(12) 


0.21939(68) 


0.33748(62) 


195 


0.22613(11) 


0.18817(69) 


0.29434(58) 


205 


0.19882(10) 


0.16401(68) 


0.25841(62) 


215 


0.175524(94) 


0.14363(70) 


0.22748(62) 


225 


0.155052(87) 


0.12328(63) 


0.19853(56) 


235 


0.137434(81) 


0.10832(51) 


0.17566(46) 


245 


0.121948(76) 


0.09545(49) 


0.15595(46) 



'/jet 




^NLO.cxcl rn b -| 




-4.32 


0.15890(34) 


0.2762(36) 


0.3316(37) 


-3.96 


0.50253(69) 


0.7129(77) 


0.8899(75) 


-3.6 


1.1812(11) 


1.463(16) 


1.891(13) 


-3.24 


2.3948(16) 


2.654(15) 


3.574(14) 


-2.88 


4.3033(22) 


4.382(22) 


6.085(22) 


-2.52 


6.9684(29) 


6.676(29) 


9.489(26) 


-2.16 


10.3269(35) 


9.476(34) 


13.772(32) 


-1.8 


14.2214(41) 


12.626(51) 


18.683(44) 


-1.44 


18.2536(47) 


16.030(49) 


23.865(38) 


-1.08 


21.9552(52) 


19.196(57) 


28.692(41) 


-0.72 


24.9131(55) 


21.607(62) 


32.554(46) 


-0.36 


26.8378(57) 


23.234(60) 


34.975(43) 





27.4891(58) 


23.82(11) 


35.877(60) 


0.36 


26.8378(57) 


23.238(63) 


35.063(54) 


0.72 


24.9131(55) 


21.608(62) 


32.531(49) 


1.08 


21.9552(52) 


19.171(62) 


28.672(39) 


1.44 


18.2536(47) 


16.002(54) 


23.846(46) 


1.8 


14.2214(41) 


12.664(52) 


18.653(51) 


2.16 


10.3269(35) 


9.480(36) 


13.773(32) 


2.52 


6.9684(29) 


6.681(29) 


9.502(23) 


2.88 


4.3033(22) 


4.386(22) 


6.085(18) 


3.24 


2.3948(16) 


2.649(14) 


3.575(13) 


3.6 


1.1812(11) 


1.452(14) 


1.888(13) 


3.96 


0.50253(69) 


0.7123(76) 


0.8883(67) 


4.32 


0.15890(34) 


0.2760(35) 


0.3309(35) 



ID + - 


do\a 

dfe+f!- 


nhl 
puj 


^o"nlo,cxc1 r„un 


rfo"NLO.incl r „un 
d<P e+ ,- LPDJ 


0.062831855 


21.1649 


(87) 


17.523(65) 


28.681(51) 


0.18849556 


23.3692 


(90) 


19.185(66) 


31.436(52) 


0.31415927 


24.5298 


(93) 


20.322(80) 


33.167(58) 


0.43982297 


24.6116 


(92) 


20.398(67) 


33.105(55) 


0.56548668 


24.7017 


(94) 


20.686(72) 


33.341(59) 


0.69115039 


24.8908 


(93) 


21.067(83) 


33.678(58) 


0.81681409 


25.1940 


(94) 


21.478(87) 


34.057(70) 


0.94247777 


25.5753 


(95) 


21.942(79) 


34.515(67) 


1.0681415 


26.0338 


(97) 


22.562(80) 


35.248(80) 


1.1938053 


26.6636 


(98) 


23.28(12) 


36.045(66) 


1.319469 


27.3938 


(99) 


24.404(96) 


37.224(61) 


1.4451327 


28.293 


(10) 


25.319(97) 


38.249(66) 


1.5707963 


29.363 


(10) 


26.56(13) 


39.692(82) 


1.69646 


30.595 


(11) 


27.62(10) 


41.299(74) 


1.8221238 


32.037 


(11) 


29.20(13) 


43.273(86) 


1.9477875 


33.644 


(11) 


30.99(14) 


45.258(81) 


2.0734512 


35.538 


(11) 


32.51(13) 


47.443(84) 


2.1991148 


37.658 


(12) 


34.63(14) 


50.073(96) 


2.3247785 


39.973 


(12) 


36.82(13) 


52.736(97) 


2.4504423 


42.582 


(12) 


39.30(14) 


55.74(10) 


2.5761059 


45.274 


(13) 


41.53(15) 


58.68(11) 


2.7017697 


47.992 


(13) 


43.90(15) 


61.631(96) 


2.8274333 


50.626 


(13) 


46.54(14) 


64.58(13) 


2.9530971 


52.935 


(14) 


48.63(15) 


66.89(11) 


3.0787608 


54.506 


(14) 


49.86(22) 


68.44(11) 



prw f) + - 


d(T LO | 


nhl 


^""NLO.excl r„UT 


f'°"NLO,incl | 


nhl 

JJUJ 


rfcosft,+„- ' 


</cos0 e +„- L P DJ 




-0.96 


49.745 


(16) 


48.56(13) 


66.057 


;93) 


-0.88 


43.269 


(15) 


41.75(16) 


58.03 


;i4) 


-0.8 


40.279 


(14) 


38.59(14) 


54.24 


;io) 


-0.72 


38.446 


(14) 


36.85(23) 


51.78 


;i4) 


-0.64 


37.272 


(13) 


35.21(20) 


50.40 


;n) 


-0.56 


36.580 


(13) 


34.52(16) 


49.37 


;i4) 


-0.48 


36.217 


(13) 


33.49(20) 


48.74 


;i4) 


-0.4 


36.137 


(13) 


33.93(22) 


48.85 


;i3) 


-0.32 


36.306 


(13) 


33.55(17) 


48.85 


;i2) 


-0.24 


36.707 


(13) 


33.78(17) 


49.43 


;i3) 


-0.16 


37.307 


(14) 


33.87(19) 


50.07 


;i2) 


-0.08 


38.142 


(14) 


34.87(23) 


51.31 


;i5) 





39.220 


(14) 


35.49(22) 


52.32 


;i3) 


0.08 


40.626 


(14) 


36.29(24) 


53.91 


;i5) 


0.16 


42.353 


(15) 


37.99(16) 


56.45 


;i2) 


0.24 


44.441 


(15) 


39.50(18) 


58.95 


;i2) 


0.32 


47.056 


(16) 


42.28(28) 


62.46 


;i3) 


0.4 


50.318 


(16) 


44.80(18) 


66.58 


;n) 


0.48 


54.494 


(17) 


47.72(18) 


71.83 


;i2) 


0.56 


59.959 


(18) 


52.68(22) 


78.87 


;i8) 


0.64 


67.158 


(20) 


58.82(19) 


88.13 


;n) 


0.72 


77.089 


(22) 


67.12(18) 


100.91 


;i4) 


0.8 


91.694 


(24) 


78.82(25) 


118.88 


;is) 


0.88 


112.407 


(28) 


95.56(19) 


145.38 


;i4) 


0.96 


118.601 


(28) 


96.50(17) 


153.08 


;i7) 



£>T,miss[GeV] 

27.5 
32.5 
37.5 
42.5 
47.5 
52.5 
57.5 
62.5 
67.5 
72.5 
77.5 
82.5 
87.5 
92.5 
97.5 
102.5 
107.5 
112.5 
117.5 
122.5 
127.5 
132.5 
137.5 
142.5 
147.5 



dcr LO 
dpr,mws 



Pb 
GeV 



0.22486(12) 
0.22916(12) 
0.22730(12) 
0.21985(11) 
0.20813(11) 
0.19202(10) 
0.172431(93) 
0.150144(84) 
0.126317(75) 
0.102584(65) 
0.081294(56) 
0.062423(47) 
0.045999(38) 
0.032865(30) 
0.023546(23) 
0.017363(19) 
0.013125(16) 
0.010179(14) 
0.007999(12) 
0.006360(10) 
0.0051262(87) 
0.0041637(78) 
0.0034107(68) 
0.0028220(62) 
0.0023492(55) 



<^°"NLO 
£*PT,mii 



Pb 
GeV 



0.26071(41) 
0.26346(38) 
0.25967(42) 
0.25002(38) 
0.23658(42) 
0.21661(33) 
0.19402(36) 
0.16747(28) 
0.13992(24) 
0.11164(24) 
0.08633(21) 
0.06500(21) 
0.04640(14) 
0.031666(87) 
0.021567(68) 
0.014919(50) 
0.010781(42) 
0.007943(40) 
0.006024(30) 
0.004588(26) 
0.003589(36) 
0.002807(18) 
0.002229(17) 
0.001784(15) 
0.001410(14) 



.ir 



Pb 

GeV 



0.29070(40) 
0.29475(38) 
0.29130(41) 
0.28133(36) 
0.26677(42) 
0.24525(33) 
0.22024(35) 
0.19098(28) 
0.16047(25) 
0.12928(25) 
0.10102(20) 
0.07693(19) 
0.05614(13) 
0.039498(86) 
0.027789(68) 
0.019946(52) 
0.014759(40) 
0.011207(34) 
0.008648(30) 
0.006727(26) 
0.005334(35) 
0.004260(17) 
0.003421(17) 
0.002783(14) 
0.002263(13) 



y i jet l * j 


dcr L0 
rf PT,jct 


r P b " 

GeV 




d<T NLO.excl 
rf PT,jct 


r P b " 

GeV 




dcr NLO.incl 
rf PT,jct 


r P b " 

GeV 




25 


0.38977(12) 


0.49226(78) 


0.51143(76) 


35 


0.214392(77) 


0.25162(56) 


0.28034(60) 


45 


0.133866(56) 


0.14689(38) 


0.17183(41) 


55 


0.089710(43) 


0.09317(33) 


0.11298(29) 


65 


0.063002(34) 


0.06205(22) 


0.07802(20) 


75 


0.045775(28) 


0.04275(21) 


0.05562(18) 


85 


0.034102(23) 


0.03022(16) 


0.04059(14) 


95 


0.025834(19) 


0.02171(20) 


0.02992(13) 


105 


0.019880(16) 


0.01599(11) 


0.02268(13) 


115 


0.015488(14) 


0.01182(10) 


0.01727(11) 


125 


0.012117(12) 


0.009083(79) 


0.013355(82) 


135 


0.009580(10) 


0.006684(69) 


0.010184(72) 


145 


0.0075986(91) 


0.005089(91) 


0.007916(67) 


155 


0.0060508(79) 


0.004025(52) 


0.006279(38) 


165 


0.0048493(69) 


0.002966(56) 


0.004837(49) 


175 


0.0038996(61) 


0.002313(35) 


0.003824(34) 


185 


0.0031335(53) 


0.001824(38) 


0.003053(36) 


195 


0.0025325(47) 


0.001419(38) 


0.002369(30) 


205 


0.0020389(41) 


0.000985(62) 


0.001861(43) 


215 


0.0016544(36) 


0.000903(30) 


0.001537(37) 


225 


0.0013444(32) 


0.000639(17) 


0.001148(17) 


235 


0.0010868(28) 


0.000527(13) 


0.000950(14) 


245 


0.0008872(25) 


0.000398(14) 


0.000729(14) 


255 


0.0007201(21) 


0.000306(18) 


0.000580(19) 


265 


0.0005831(19) 


0.000256(13) 


0.000473(13) 



[GeV] 

27.5 
32.5 
37.5 
42.5 
47.5 
52.5 
57.5 
62.5 
67.5 
72.5 
77.5 
82.5 
87.5 
92.5 
97.5 
102.5 
107.5 
112.5 
117.5 
122.5 
127.5 
132.5 
137.5 
142.5 
147.5 



Pb 

GeV 



0.30973(16) 
0.30491(14) 
0.27515(13) 
0.23554(11) 
0.195422(96) 
0.159132(84) 
0.128552(74) 
0.103607(65) 
0.083574(58) 
0.067426(51) 
0.054683(45) 
0.044524(40) 
0.036523(36) 
0.029974(32) 
0.024674(29) 
0.020515(26) 
0.017003(23) 
0.014195(21) 
0.011881(19) 
0.009974(17) 
0.008367(16) 
0.007092(14) 
0.006007(13) 
0.005081(12) 
0.004332(11) 



do nlo 



Pb 
GeV 



do~ NLOjnc 



Pb 

GeV 



0.35406(48) 
0.35174(44) 
0.31884(49) 
0.26988(49) 
0.22188(45) 
0.17840(26) 
0.14200(28) 
0.11224(22) 
0.08928(25) 
0.07061 (24) 
0.05642(16) 
0.04495(15) 
0.03607(23) 
0.02967(22) 
0.023881(94) 
0.019219(88) 
0.015744(86) 
0.012863(68) 
0.010706(63) 
0.008755(49) 
0.007244(48) 
0.006093(39) 
0.004931(36) 
0.004208(33) 
0.003622(38) 



0.39974(48) 
0.39591 (45) 
0.35935(49) 
0.30501 (47) 
0.25183(43) 
0.20373(26) 
0.16308(26) 
0.12977(23) 
0.10379(27) 
0.08269(25) 
0.06644(16) 
0.05324(13) 
0.04300(20) 
0.03551(21) 
0.028747(94) 
0.023307(87) 
0.019181(86) 
0.015729(67) 
0.013152(64) 
0.010832(52) 
0.009042(53) 
0.007605(39) 
0.006200(36) 
0.005308(36) 
0.004554(39) 



/?T,e+[GeV] 

27.5 
32.5 
37.5 
42.5 
47.5 
52.5 
57.5 
62.5 
67.5 
72.5 
77.5 
82.5 
87.5 
92.5 
97.5 
102.5 
107.5 
112.5 
117.5 
122.5 
127.5 
132.5 
137.5 
142.5 
147.5 



dp Ts + 



Pb 

GeV 



0.30969(16) 
0.30476(14) 
0.27514(12) 
0.23547(11) 
0.195342(96) 
0.159304(85) 
0.128438(74) 
0.103669(66) 
0.083655(58) 
0.067429(51) 
0.054722(46) 
0.044628(41) 
0.036517(36) 
0.029963(33) 
0.024763(29) 
0.020480(26) 
0.017023(24) 
0.014164(21) 
0.011861(19) 
0.009959(17) 
0.008412(16) 
0.007077(14) 
0.006009(13) 
0.005095(12) 
0.004309(11) 



excj 
rf PT,c+ 



Pb 

GeV 



do NLO.incl 



Pb 
GeV 



0.35540(53) 
0.35176(47) 
0.31748(41) 
0.27055(38) 
0.22209(32) 
0.17864(32) 
0.14141(25) 
0.11270(25) 
0.08918(20) 
0.07026(24) 
0.05676(21) 
0.04505(14) 
0.03602(11) 
0.02940(10) 
0.02358(10) 
0.019326(99) 
0.015724(80) 
0.012928(72) 
0.010682(67) 
0.008731(52) 
0.007133(66) 
0.006146(67) 
0.005094(46) 
0.004191(34) 
0.003517(34) 



0.40097(53) 
0.39624(46) 
0.35780(41) 
0.30585(37) 
0.25213(33) 
0.20402(31) 
0.16253(25) 
0.13015(23) 
0.10373(19) 
0.08238(26) 
0.06667(20) 
0.05339(14) 
0.04290(11) 
0.03518(10) 
0.028413(98) 
0.023401 (98) 
0.019199(83) 
0.015834(66) 
0.013163(68) 
0.010833(52) 
0.008860(68) 
0.007656(66) 
0.006365(45) 
0.005288(33) 
0.004450(34) 



daro 



[pb] 



,excl 



[pb] 



^Q"NLO.ir 



[pb] 



-2.4 
-2.2 
-2 
-1.8 
-1.6 
-1.4 
-1.2 
-1 
-0.8 
-0.6 
-0.4 
-0.2 

0.2 
0.4 
0.6 
0.8 
1 

1.2 
1.4 
1.6 
1.8 
2 
2.2 
2.4 



0.27168(42) 
0.46828(63) 
0.74839(87) 
1.1136(11) 
1.5503(15) 
2.0437(17) 
2.5502(20) 
3.0353(22) 
3.4528(24) 
3.7761(25) 
3.9845(25) 
4.0658(25) 
4.0259(24) 
3.8716(23) 
3.6187(21) 
3.2972(19) 
2.9122(17) 
2.5006(15) 
2.0766(13) 
1.6645(12) 
1.2817(10) 
0.94199(84) 
0.66057(69) 
0.43742(55) 
0.27079(42) 



0.3040(16) 
0.5228(60) 
0.8256(40) 
1.2392(57) 
1.6979(98) 
2.2547(85) 
2.8151(98) 
3.3262(83) 
3.782(13) 
4.130(12) 
4.345(12) 
4.442(12) 
4.4000(97) 
4.244(12) 
3.9915(95) 
3.6443(82) 
3.2341(71) 
2.7948(62) 
2.3340(57) 
1.8748(68) 
1.4585(44) 
1.0699(34) 
0.7577(31) 
0.5046(26) 
0.3175(21) 



0.3409(16) 
0.5895(62) 
0.9343(40) 
1.4061(57) 
1.948(10) 
2.5763(89) 
3.218(10) 
3.8198(86) 
4.367(12) 
4.768(10) 
5.019(12) 
5.135(12) 
5.0760(97) 
4.903(12) 
4.583(10) 
4.1792(86) 
3.6921 (68) 
3.1792(64) 
2.6390(58) 
2.1021(68) 
1.6349(48) 
1.1884(33) 
0.8373(31) 
0.5544(27) 
0.3468(21) 





L P UJ 


ATNLO,exclr 

rf^ e + L P UJ 


L P UJ 


-2.4 


0.27124(43) 


0.3139(20) 


0.3443(19) 


-2.2 


0.43757(56) 


0.5046(25) 


0.5550(25) 


-2 


0.66082(71) 


0.7610(37) 


0.8416(33) 


-1.8 


0.94463(87) 


1.0770(39) 


1.1947(39) 


-1.6 


1.2805(10) 


1.4555(44) 


1.6258(44) 


-1.4 


1.6640(12) 


1.8749(53) 


2.1066(53) 


-1.2 


2.0764(14) 


2.3277(59) 


2.6312(62) 


-1 


2.4991(16) 


2.7865(63) 


3.1675(63) 


-0.8 


2.9163(17) 


3.2388(74) 


3.6926(85) 


-0.6 


3.2962(19) 


3.6422(73) 


4.1724(88) 


-0.4 


3.6215(21) 


3.9869(91) 


4.5880(85) 


-0.2 


3.8654(23) 


4.2419(91) 


4.8908(88) 





4.0231(24) 


4.397(11) 


5.082(11) 


0.2 


4.0651(25) 


4.4434(96) 


5.1421(99) 


0.4 


3.9837(25) 


4.335(10) 


5.025(10) 


0.6 


3.7760(25) 


4.133(15) 


4.776(15) 


0.8 


3.4522(24) 


3.788(10) 


4.364(10) 


1 


3.0365(22) 


3.3306(98) 


3.8285(93) 


1.2 


2.5526(20) 


2.812(14) 


3.213(12) 


1.4 


2.0468(17) 


2.2697(98) 


2.5971(95) 


1.6 


1.5535(14) 


1.7213(84) 


1.9599(84) 


1.8 


1.1129(11) 


1.2313(55) 


1.3990(55) 


2 


0.74748(84) 


0.8308(60) 


0.9351(55) 


2.2 


0.46781(60) 


0.5242(31) 


0.5901(30) 


2.4 


0.27187(41) 


0.3043(17) 


0.3419(17) 





do~LO 

dm c +^- 


r P b " 

GeV 




dcr NLO.excl 

<2m e+/J - 


r P b " 

GeV 




dcr NLO.incl 


r P b " 

GeV 




5 


0.002587(10) 


0.003042(64) 


0.003389(59) 


15 


0.016720(24) 


0.018743(67) 


0.021646(66) 


25 


0.028051(32) 


0.031663(92) 


0.036463(91) 


35 


0.039593(38) 


0.04497(12) 


0.05156(12) 


45 


0.052892(45) 


0.06066(14) 


0.06926(14) 


55 


0.073373(54) 


0.08432(17) 


0.09572(17) 


65 


0.093311(59) 


0.10709(21) 


0.12126(21) 


75 


0.100353(59) 


0.11545 (21) 


0.13094(21) 


85 


0.096957(56) 


0.11012(20) 


0.12545(20) 


95 


0.087873(51) 


0.09854(19) 


0.11277(19) 


105 


0.076902(46) 


0.08558(16) 


0.09817(16) 


115 


0.065851(41) 


0.07263(16) 


0.08357(16) 


125 


0.055874(37) 


0.06115(13) 


0.07047(13) 


135 


0.047015(33) 


0.05058(14) 


0.05842(14) 


145 


0.039484(30) 


0.04268(16) 


0.04919(15) 


155 


0.033171(27) 


0.03519(12) 


0.04060(12) 


165 


0.027827(25) 


0.029571(97) 


0.034151(97) 


175 


0.023336(23) 


0.024458(88) 


0.028257(87) 


185 


0.019641(21) 


0.02064(11) 


0.02381(11) 


195 


0.016627(19) 


0.017056(93) 


0.019698(93) 


205 


0.013979(18) 


0.014438(92) 


0.016646(91) 


215 


0.011843(16) 


0.011997(72) 


0.013844(69) 


225 


0.010019(15) 


0.010061(62) 


0.011623(63) 


235 


0.008489(14) 


0.008514(58) 


0.009808(58) 


245 


0.007250(13) 


0.007370(52) 


0.008465(52) 



'/jet 




^0"NLO,cxcl r n l,T 


^CTNLO.incl r Li 


-2.4 


0.57315(68) 


0.6833(78) 


0.7509(86) 


-2.2 


0.78791(83) 


0.901(11) 


1.005(11) 


-2 


1.0390(10) 


1.221(16) 


1.346(16) 


-1.8 


1.3276(12) 


1.533(19) 


1.711(13) 


-1.6 


1.6330(13) 


1.866(18) 


2.083(13) 


-1.4 


1.9453(15) 


2.166(14) 


2.483(15) 


-1.2 


2.2666(17) 


2.523(13) 


2.875(14) 


-1 


2.5753(18) 


2.815(15) 


3.218(17) 


-0.8 


2.8418(19) 


3.072(15) 


3.572(16) 


-0.6 


3.0669(20) 


3.326(20) 


3.838(24) 


-0.4 


3.2377(21) 


3.488(25) 


4.097(32) 


-0.2 


3.3455(21) 


3.592(25) 


4.174(19) 





3.3903(22) 


3.630(28) 


4.219(21) 


0.2 


3.3544(22) 


3.586(21) 


4.191(18) 


0.4 


3.2405(22) 


3.482(18) 


4.070(19) 


0.6 


3.0610(21) 


3.307(17) 


3.829(18) 


0.8 


2.8506(21) 


3.111(20) 


3.582(21) 


1 


2.5775(19) 


2.816(17) 


3.240(17) 


1.2 


2.2677(18) 


2.480(19) 


2.841(17) 


1.4 


1.9445(16) 


2.180(16) 


2.482(16) 


1.6 


1.6339(15) 


1.833(13) 


2.086(13) 


1.8 


1.3294(13) 


1.516(12) 


1.700(13) 


2 


1.0412(11) 


1.195(13) 


1.338(11) 


2.2 


0.78676(95) 


0.9427(94) 


1.0262(85) 


2.4 


0.57316(77) 


0.6867(89) 


0.7516(84) 



fe + fi- 

0.062831855 
0.18849556 
0.31415927 
0.43982297 
0.56548668 
0.69115039 
0.81681409 
0.94247777 
1.0681415 
1.1938053 
1.319469 
1.4451327 
1.5707963 
1.69646 
1.8221238 
1.9477875 
2.0734512 
2.1991148 
2.3247785 
2.4504423 
2.5761059 
2.7017697 
2.8274333 
2.9530971 
3.0787608 



dcr L0 
dip e+ - 



[pb] 



1.5811(19) 
1.7229(20) 
1.8202(20) 
1.8724(21) 
1.9321(21) 
2.0126(21) 
2.1161(22) 
2.2367(22) 
2.3776(23) 
2.5414(23) 
2.7318(24) 
2.9411(26) 
3.1711(26) 
3.4213(27) 
3.6860(27) 
3.9508(28) 
4.2263(29) 
4.5080(30) 
4.7721(31) 
5.0090(32) 
5.2482(33) 
5.4610(34) 
5.6843(34) 
5.9002(35) 
6.0086(36) 



,exi 

1.7284 
1.8614 
1.9833 
2.036 
2.1060 
2.2103 
2.3300 
2.4762 
2.6276 
2.8411 
3.0415 
3.299 
3.584 
3.850 
4.135 
4.436 
4.722 
5.053 
5.309 
5.540 
5.821 
5.983 
6.205 
6.390 
6.479 



[ [pb] 

(86) 
(80) 
(98) 
(14) 
(73) 
(85) 
(99) 
(94) 
(98) 
(95) 
(96) 

(11) 
(20) 

(11) 
(12) 
(12) 
(12) 
(15) 
(12) 
(12) 
(14) 
(13) 
(12) 
(13) 
(12) 



dcr NLO ,incl r —1-1 

2.0408(86) 
2.2075(80) 
2.340(10) 
2.405(13) 
2.4815(72) 
2.5944(84) 
2.719(10) 
2.8831(90) 
3.048(10) 
3.2849(99) 
3.5011(98) 
3.782(11) 
4.082(18) 
4.394(11) 
4.710(11) 
5.052(12) 
5.368(11) 
5.727(15) 
6.023(12) 
6.288(13) 
6.590(14) 
6.782(13) 
7.041(12) 
7.240(12) 
7.339(12) 



cos 9. 



-0.96 
-0.88 

-0.8 
-0.72 
-0.64 
-0.56 
-0.48 

-0.4 
-0.32 
-0.24 
-0.16 
-0.08 


0.08 
0.16 
0.24 
0.32 

0.4 
0.48 
0.56 
0.64 
0.72 

0.8 
0.88 
0.96 



rfcos ft, 



T fl 



8.9561(59) 
8.0391(55) 
7.4532(53) 
7.0061(50) 
6.6385(49) 
6.3402(47) 
6.0737(46) 
5.8496(45) 
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5.4969(44) 
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5.2193(42) 
5.0918(41) 
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4.7384(38) 
4.6672(37) 
4.5956(36) 
4.5262(36) 
4.4354(35) 
4.3315(34) 
4.1823(32) 
3.9323(31) 
3.2732(27) 
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dcosft, +„- LP°J 



9.712(22) 
8.789(26) 
8.214(24) 
7.726(29) 
7.319(24) 
7.054(26) 
6.717(33) 
6.478(33) 
6.269(21) 
6.139(31) 
5.932(37) 
5.801(29) 
5.647(20) 
5.577(34) 
5.447(23) 
5.350(30) 
5.249(21) 
5.175(18) 
5.078(22) 
5.003(18) 
4.905(18) 
4.780(20) 
4.581(16) 
4.308(13) 
3.539(12) 



10.986(21) 
9.964(23) 
9.309(27) 
8.794(26) 
8.324(23) 
8.024(29) 
7.674(30) 
7.353(32) 
7.179(22) 
7.015(29) 
6.787(31) 
6.658(26) 
6.489(25) 
6.387(34) 
6.256(22) 
6.168(30) 
6.041(21) 
5.955(18) 
5.847(21) 
5.764(19) 
5.663(19) 
5.534(23) 
5.311(16) 
5.018(13) 
4.151(12) 



[9] S. Dittmaier, S. Kallweit and P. Uwer, Phys. Rev. Lett. 100 (2008) 062003 
HarXiv:0710.1577 [hep-ph]]. 

[10] J. M. Campbell, R. Keith Ellis and G. Zanderighi, JHEP 0712 (2007) 056 
HarXiv:0710.1832l [hep-ph]] . 



[11] G. Sanguinetti and S. Karg. arXiv:0806. 13941 rhep-phl ; 

T. Binoth, J. -P. Guillet, S. Karg, N. Kauer and G. Sanguinetti, in preparation. 

[12] T. Kobayashi, Prog. Theor. Phys. 49 (1973) 282. 

[13] N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531. 

[14] C. Neu [CDF Collaboration], FERMILAB-CONF-06-162-E. 

[15] K. Hanagaki [DO Collaboration], FERMILAB-CONF-05-647-E. 

[16] T. Wright [CDF and DO Collaboration! . larXiv:0707J7T2l fhep-exl . 

[17] A. Wildauer and F. Akesson, in "Interlaken 2004, Computing in high energy physics and 
nuclear physics," 419-422. 

[18] M. Lehmacher. larXiv:0809.4896l fhep-exl . 

[19] A. Denner, S. Dittmaier, M. Roth and D. Wackeroth, Nucl. Phys. B 560 (1999) 33 
[ [arXiv:hep-ph/9904472l . 



[20] A. Denner, S. Dittmaier, M. Roth and L. H. Wieders, Nucl. Phys. B 724 (2005) 247 
[ [arXiv:hep-ph/0505042l . 



[21] A. Brandenburg and P. Uwer, Nucl. Phys. B 515, 279 (1998) [arXiv:hep-ph/9708350| 



[22] A. Banfi, G. P. Salam and G. Zanderighi, Eur. Phys. J. C 47 (2006) 113 



[arXiv:hep-ph/0601139| 



[23] W. Beenakker, S. Dittmaier, M. Kramer, B. Plumper, M. Spira and P. M. Zerwas, Nucl. 
Phys. B 653 (2003) 151 [arXiv:hep-ph/0211352| . 

[24] S. Dittmaier, P. Uwer and S. Weinzierl, Phys. Rev. Lett. 98 (2007) 262002 
[|arXiv:hep-ph/0703120] and Eur. Phys. J. C 59 (2009) 625 HarXiv:08 10.04521 rhep-phl 1 . 



[25] J. Kiiblbeck, M. Bohm and A. Denner, Comput. Phys. Commun. 60 (1990) 165. 



[26] S. Dittmaier, Nucl. Phys. B 675 (2003) 447 [arXiv:hep-ph/0308246| 



[27] A. Denner and S. Dittmaier, Nucl. Phys. B 658 (2003) 175 [arXiv:hep-ph/02 12259] 



[28] G. Passarino and M. J. G. Veltman, Nucl. Phys. B 160 (1979) 151. 
[29] G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B 153, 365 (1979). 
[30] W. Beenakker and A. Denner, Nucl. Phys. B 338 (1990) 349. 

54 



[31 

[32 
[33 

[34 



[35 
[36 

[37 
[38 
[39 

[40 

[41 
[42 

[43 
[44 

[45 
[46 
[47 
[48 

[49 
[50 



A. Denner, U. Nierste and R. Scharf, Nucl. Phys. B 367, 637 (1991). 

A. Denner and S. Dittmaier, Nucl. Phys. B 734, 62 (2006) [ |arXiv:hep-ph/0509 1411 . 



T. Binoth, J. P. Guillet, G. Heinrich, E. Pilon and C. Schubert, JHEP 0510 (2005) 015 
[ [arXiv:hep-ph/05042671 . 



W. Giele, E. W. N. Glover and G. Zanderighi, Nucl. Phys. Proc. Suppl. 135 (2004) 275 



[arXiv:hep-ph/0407016|; 



R. K. Ellis, W. T. Giele and G. Zanderighi, Phys. Rev. D 73 (2006) 014027 



[arXiv:hep-ph/0508308|. 



T. Hahn, Comput. Phys. Commun. 140 (2001) 418 [ |arXiv:hep-ph/00 1 22601 • 

T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118 (1999) 153 
]arXiv:hep-ph/9807565| . 

G. J. van Oldenborgh and J. A. M. Vermaseren, Z. Phys. C 46, 425 (1990). 
G. J. van Oldenborgh, Comput. Phys. Commun. 66 (1991) 

Z. Bern, L. J. Dixon and D. A. Kosower, Nucl. Phys. B 412 (1994) 751 
garXiv:hep-ph/93062401 . 

G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B 44 (1972) 189; 

P. Breitenlohner and D. Maison, Commun. Math. Phys. 52 (1977) 39. 



F. Jegerlehner, Eur. Phys. J. C 18 (2001) 673 [arXiv:hep-th/0005255| 



S. Catani and M. H. Seymour, Nucl. Phys. B 485 (1997) 291 [Erratum-ibid. B 510 (1998) 
503] [arXiv:hep-ph/9605323| . 



S. Dittmaier, Phys. Rev. D 59 (1999) 016007 [ |arXiv:hep- ph/98054451. 

F. A. Berends, R. Pittau and R. Kleiss, Nucl. Phys. B 424 (1994) 308 
|arXiv:hep-ph/9404313|. 



R. Kleiss and R. Pittau, Comput. Phys. Commun. 83, 141 (1994) [arXiv:hep-ph/9405257| 



M. Roth, doctoral thesis, DISS-ETH-13363, 1999, arXiv:hep-ph/0008033 



W. Kilian, T. Ohl and J. Reuter. larXiv:0708.4233 [hep-ph]. 

T. Gleisberg, S. Hoche, F. Krauss, A. Schalicke, S. Schumann and J. C. Winter, JHEP 0402 
(2004) 056 [arXiv:hep-ph/0311263| . 

S. Kallweit, diploma thesis (in German), LMU Munich, 2006. 

T. Stelzer and W. F. Long, Comput. Phys. Commun. 81 (1994) 357 



[arXiv:hep-ph/9401258|. 



[51] S. Eidelman et al. [Particle Data Group], Phys. Lett. B 592 (2004) 1. 

55 



[52] A. Djouadi, J. Kalinowski and M. Spira, Comput. Phys. Commun. 108 (1998) 56 
| )hep-ph/9704448) . 

[53] S. Dittmaier and M. Kramer, Phys. Rev. D 65 (2002) 073007 [ |arXiv:hep-ph/0 1 09062[ | . 

[54] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667 (2008) 1. 

[55] S. Dittmaier, S. Kallweit and P. Uwer, PoS RADCOR2007 (2007) 009 HarXiv:0712.3163l 
[hep-ph]]. 

[56] S. Dittmaier, S. Kallweit and P. Uwer, Nucl. Phys. Proc. Suppl. 183 (2008) 279 
HarXiv:0807.0708l [hep-ph] ] . 

[57] R. K. Ellis, W. J. Stirling, B. R. Webber, "QCD and Collider Physics," Cambridge Uni- 
versity Press, 1996. 

[58] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207 
(2002) 012 garXiv:hep-ph/02011951 . 

[59] D. Stump, J. Huston, J. Pumplin, W. K. Tung, H. L. Lai, S. Kuhlmann and J. F. Owens, 
JHEP 0310 (2003) 046 [ |arXiv:hep-ph/0303013 |. 

[60] A. D. Martin, W. J. Stirling and R. S. Thorne, Phys. Lett. B 636 (2006) 259 
| |arXiv:hep-ph/0603 143| . 

[61] S. D. Ellis and D. E. Soper, Phys. Rev. D 48 (1993) 3160 [arXiv:hep-ph/9305266| . 

[62] L. J. Dixon, Z. Kunszt and A. Signer, Phys. Rev. D 60 (1999) 114037 
HarXiv:hep-ph/9907305| . 

[63] M. Dittmar and H. K. Dreiner, |arXiv:hep-ph/970340 1 [ 



56 



